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PREFACE 


The  work  described  herein  was  performed  by  Dr.  Ronald  A.  Erchul  during 
the  period  from  June  1986  to  July  1987  under  a  contract  with  the  Virginia 
Military  Institute  (VMI) ,  Lexington,  VA.  The  work  was  authorized  by  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  MS,  under 
Contract  No.  DACW39-86-K-0006 .  The  work  was  funded  by  the  Repair,  Evaluation 
Maintenance,  and  Rehabilitation  (REMR)  Research  Project  "Geophysical  Tech¬ 
niques  for  Assessment  of  Existing  Structural  Foundations,"  Work  Unit 
No.  WU  1-32315.  This  is  the  second  of  two  reports  in  a  series  of  three. 

The  field  work  was  performed  during  the  period  from  July  1986  to  May 
1987  by  Messrs.  Dennis  W.  Slifer,  Robert  Thren,  Steven  Richards,  and 
Cadets  D.  R.  Gibblngs  and  R.  W.  Marsh,  VMI.  The  data  analysis  phase  of  the 
investigation  was  performed  by  Mr.  Dennis  W.  Slifer  and  Cadets  D.  R.  Gibbings 
and  R.  W.  Marsh.  This  report  was  written  by  Dr.  Robert  A.  Erchul  and 
Mr.  Dennis  W.  Slifer.  Principal  Investigators  for  the  REMR  research  project 
are  Dr,  Dwain  K.  Butler  and  Mr.  Jose  L.  Llopis,  Earthquake  Engineering  and 
Geophysics  Division,  Geotechnical  Laboratory,  WES.  Mr.  Benjamin  I.  Kelley, 
Headquarters,  US  Army  Corps  of  Engineers,  was  Technical  Monitor  for  this  work 
The  report  was  edited  by  Ms.  Joyce  H.  Walker,  Information  Products  Division, 
Information  Technology  Laboratory,  WES. 

Acting  Commander  and  Director  of  WES  during  preparation  of  this 
report  was  LTC  Jack  R.  Stephens,  EN.  Dr.  Robert  W.  Whalin  is  Technical 
Director . 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (met¬ 
ric)  units  as  follows: 


_ Multiply _ 

degrees  (angle) 

Fahrenheit  degrees 

feet 

gallons 

inches 

miles  (US  Statute) 
milllmbus  per  foot 
ohm-foot 

pounds  (force)  seconds 
per  square  foot 

pounds  (force)  per 
square  foot 


.  By _ 

0.01745329 

5/9 

0.3048 

3.785412 

25.4 

1.609347 

3.28 

0.3048 

14.88164 

47.88026 


_ To  Obtain _ 

radians 

Celsius  degrees  or  Kelvins* 
metres 

cubic  decimetres 

millimetres 

kilometres 

milllmbus  per  metre 

ohm-metres 

poise 

pascals 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  =  (5/9)  (F  -  32).  To  obtain  Kelvin  (K)  read¬ 
ings,  use:  K  *  (5/9)  (F  -  32)  +  273.15. 
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GEOTECHNICAL  APPLICATIONS  OF  THE  SELF  POTENTIAL  (SP)  METHOD 


The  Use  of  Self  Potential  to  Detect  Groundwater  Flow  in  Karst 


PART  I:  INTRODUCTION 

Background 

1.  Report  1  in  this  series  (Erchul  1987)  presented  a  preliminary  proof- 
of-concept  field  evaluation  of  a  technique  for  mapping  subsurface  water  flow 
into/from  sinkholes  and  between  sinkholes.  The  evaluation  used  a  novel  scheme 
consisting  of  (a)  circumferential  self  potential  (SP)  electrode  arrays  to 
detect  flow  into  or  from  sinkholes  and  then  of  (b)  rectangular  SP  electrode 
arrays  to  track  or  map  flow  paths  "cross  country"  between  sinkholes.  Results 
of  the  proof-of-concept  field  evaluation  were  positive  and  indicated  that  the 
concept  might  lead  to  a  simple,  inexpensive  procedure  for  mapping  subsurface 
water  flow  paths  in  karst  areas.  This  current  report  (Report  2)  extends  the 
preliminary  work  to  investigate: 

a.  Environmental  effects  (temperature,  rain,  vegetation,  etc.)  on 
SP  measurements  with  metal  electrodes. 

b.  Relationships  between  flow  quantity  and  SP  magnitudes. 

c.  Limitations  of  method. 

Review 

2.  Karst  regions,  areas  of  soluble  carbonate  bedrock,  have  unique 
hydrological  characteristics.  Environmental  and  engineering  problems  result 
from  the  solutional  removal  of  bedrock — a  process  that  produces  sinkholes, 
preferential  subsurface  drainage,  and  random  soil  depths.  Water  supply  and 
quality,  drainage,  flooding,  and  land  use  problems  are  different  from  those  in 
nonkarst  areas.  Karst  terrain  comprises  approximately  twenty  percent  of  the 
land  area  of  the  United  States  and  is  extensive  worldwide  (Quinlan  and  Ewers 
1986) .  Although  this  research  was  conducted  in  the  karst  of  the  Valley  and 
Ridge  Province  of  the  Appalachian  region,  the  results  should  be  applicable  to 
most  other  karst  regions  and  perhaps  nonkarst  areas  as  well. 
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3.  There  is  always  a  need  for  more  information  about  subsurface  condi¬ 
tions  and  processes  in  the  earth.  The  SP  technique  promises  to  be  an  effec¬ 
tive  method  for  producing  such  information.  Although  it  is  one  of  the  least 
understood  geophysical  methods,  recent  research  and  applications  indicate  that 
because  of  a  close  relationship  between  SP  and  hydrogeologic  parameters,  SP 
may  become  a  valuable  and  standard  technique  in  hydrogeology  (Ernston  and 
Scherer  1986) . 

4.  SP  is  an  electrical  geophysical  technique  that  measures  naturally 
occurring  voltage  fields  in  the  earth's  surface.  SP  is  one  of  the  oldest  geo¬ 
physical  methods,  first  used  by  Robert  Fox  in  the  copper  mines  of  Cornwall, 
England  in  1830.  As  a  prospecting  tool  for  sulphide  ores  ana  graphite  depos¬ 
its,  SP  has  been  widely  used  by  geologists  around  the  world.  Likewise,  SP 
logs  are  commonly  made  from  boreholes  and  wells  to  record  changes  in  perme¬ 
ability,  pore  fluid,  and  lithology. 

5.  In  recent  years,  the  SP  technique  has  received  Increased  study  for  a 
variety  of  novel  geotechnical  applications.  Geothermal  exploration  using  SP 
has  been  successfully  demonstrated  in  many  locations  (Anderson  and  Johnson 
1976;  Corwin  and  Hoover  1978).  SP  has  also  been  used  to  detect  and  track 
underground  coal  mine  fires  (Rodriguez  1984;  Corwin  and  Hoover  1978).  Elec¬ 
tromagnetic  anomalies  induced  by  water  diffusion  prior  to  earthquakes  have 
been  studied  by  Fltterman  (1978)  and  Ishido  (1981).  SP  monitoring  techniques 
in  this  application  may  eventually  produce  a  useful  predictive  tool.  Soviet 
scientists  (Ogilvy,  Ayed,  and  Bogoslovsky  1969;  Bogoslovsky  and  Ogilvy  1972) 
were  among  the  first  to  experiment  with  SP  methods  in  applied  hydrological  and 
geotechnical  uses.  They  found  that  SP  was  able  to  detect  leakages  at  dams, 
reservoirs,  and  drainage  structures,  and  it  is  this  application  that  has 
potential  for  REMR  applications.  Another  useful  application  is  evaluating 
ground-water  patterns  relating  to  landslide  processes. 

6.  In  the  United  States,  the  application  of  SP  to  geotechnical  problems 
has  been  studied  by  the  US  Army  Corps  of  Engineers  at  some  of  its  dam  sites. 
Cooper,  Koester  and  Franklin  (1982)  conducted  an  SP  survey  at  Gathright  Dam  in 
Virginia.  Other  studies  Include  Clearwater  Dam  in  Missouri  (Koester,  Butler, 
Cooper,  and  Llopis  1984),  Mill  Creek  Dam  in  Washington  (Butler,  Wahl,  and 
Sharp  1984),  Medford  Cave  and  Manatee  Springs  in  Florida  (Cooper  1983),  and 
Beaver  Dam  in  Arkansas  (Llopis  and  Butler  1988).  Other  agencies  (US  Bureau  of 
Reclamation)  and  consulting  firms  are  beginning  to  use  SP 
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methods  for  seepage  analysis  and  to  define  hazardous  areas  In  landslide  masses 
(Marklewlez  and  Randall  1984) . 

7.  Many  of  the  studies  conducted  by  the  Corps  of  Engineers  are  of  par¬ 
ticular  Interest  because  of  their  locations  In  karst  terrain.  Impoundments  In 
karst  areas  have  always  been  plagued  by  leakage  problems  because  of  the  pres¬ 
ence  of  caves >  sinkholes,  erratic  and  heterogeneous  ground-water  patterns,  and 
other  solutlonal  features  which  contribute  to  unpredictable  subsurface  condi¬ 
tions.  It  Is  not  necessary  to  stop  all  leakage,  but  It  Is  necessary  to  con¬ 
trol  It  to  the  extent  that  there  is  no  danger  of  piping  or  uplift-induced 
failure.  To  control  leakage  from  an  existing  water  resources  structure  or  its 
foundation,  one  has  first  to  find  the  seepage  path  or  paths. 

8.  Karst  terrain  is  especially  vulnerable  to  ground-water  contamina¬ 
tion.  Researchers  have  used  other  geophysical  techniques  (electrical  resis¬ 
tivity,  seismic  refraction,  gravimetry,  terrain  conductivity,  ground  probing 
radar,  etc.)  In  studying  subsurface  features  of  karst,  but  SF  Is  attractive  as 
a  cost-effective  method  for  Initial  investigations  of  an  area.  This  is  due  In 
part  to  the  use  of  simple  metal  electrodes  which  was  pioneered  on  Corps  of 
Engineers  projects;  this  adaptation  makes  the  technique  relatively  fast  and 
easy  to  execute.  This  simplification  of  the  technique  combined  with  the 
apparent  environmental  and  engineering  applications  In  sensitive  karst  areas 
are  the  basis  for  the  present  study.  As  a  simple,  fast,  and  Inexpensive 
method,  SP  should  be  valuable  for  augmenting  other  geophysical  and  geotech¬ 
nical  subsurface  explorative  methods.  It  seems  especially  well  suited  for 
monitoring  changes  over  relatively  long  periods  of  time.  In  addition  to  its 
utility  for  helping  to  solve  site-specific  geotechnical  problems,  SP  research 
in  karst  areas  may  contribute  to  advances  in  theoretical  concepts  about 
ground-water  movement,  storage,  and  recharge  processes. 

9.  The  purpose  of  this  research  Is  to  evaluate  the  SP  technique  (as 
developed  by  the  Corps  of  Engineers)  by  applying  it  to  a  monitoring  program  at 
a  variety  of  karst  features  within  a  study  area  In  west  central  Virginia.  The 
SP  results  are  evaluated  by  comparison  with  site  observations,  geological 
studies,  results  from  other  geophysical  methods,  stream  flow  measurements, 
cave  surveys,  artificial  water  Injections  (recharge)  into  sinkholes,  dye  trac¬ 
ing,  and  monitoring  of  soil  temperature,  precipitation,  ground-water  conduc¬ 
tivity,  and  other  environmental  parameters  that  may  affect  SP.  In  addition,  a 
computerized  data  collection  system  Is  developed  for  field  use.  This  system 
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greatly  Increases  the  capabilities  and  effectiveness  of  monitoring  both 
short-term  hydrologic  recharge  events  as  well  as  long-term  and  seasonal 
changes  at  the  study  sites. 


Theory  of  Streaming  Potential 

10.  As  used  in  geophysics,  SP  refers  to  a  phenomenon  that  originates 
from  several  different  mechanisms.  Electrochemical  SP  effects  are  caused  by 
differences  in  solution  concentrations  and  oxidation  states  of  earth  mate¬ 
rials.  This  is  the  basis  for  the  earth  cell  or  galvanic  (battery)  cell  expla¬ 
nation  commonly  used  in  connection  with  SP  anomalies  generated  by  metallic  ore 
bodies.  The  electroklnetic  (or  electrofiltration)  mechanism,  on  the  other 
hand,  generates  potentials  by  the  flow  of  water  or  heat  through  a  porous 
medium.  SP  originating  from  this  source  is  known  as  streaming  potential  and 
is  of  primary  Interest  in  ground-water  investigations  and  the  present  study. 

11.  Electrical  effects  resulting  from  the  application  of  man-made  cur¬ 
rents  to  the  earth  are  not  Included  in  the  definition  of  SP,  hence  the  term 
"spontaneous."  Telluric  currents  (elecuromagnetlcally  induced  earth  currents) 
are  also  excluded  from  a  consideration  of  SP  effects.  This  is  not  to  say  that 
cultural  or  telluric  signals  can  be  ignored  in  SP  work;  in  some  situations 
they  can  produce  noise  which  must  be  recognized  and  filtered  out.  Another 
possible  source  for  signal  noise  in  SP  work  may  derive  from  "bioelectric  cur¬ 
rents"  from  vegetation  (Nourbehecht  1963).  Presumably,  this  effect  can  be 
recognized  as  obvious  differences  in  background  between  wooded  and  open  ter¬ 
rain  when  conducting  traverses  across  such  boundaries.  Electrode  polarization 
effects  can  also  Influence  SP  results.  Electrodes  and  techniques  of  use  have 
been  debated  in  the  literature.  The  relatively  recent  trend  toward  using 
metal  electrodes  (as  opposed  to  porous  ceramic  nonpolarizing  electrodes)  is 
still  somewhat  controversial.  More  detail  on  this  topic  is  given  later  in 
this  study. 

12.  Interpretation  of  the  results  from  the  present  study  is  based  on 
the  detection  of  negative  anomalies  in  the  potential  field  which  are  assumed 
to  be  primarily  due  to  streaming  potential  effects.  Therefore,  an  increasing 
SP  value  is  one  that  is  becoming  more  negative;  whereas,  decreasing  SP  value 
is  one  that  is  becoming  less  negative  (or  more  positive).  Ogilvy  (1967) 
states  that  the  most  significant  results  are  obtained  in  karst  regions  where 
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considerable  quantities  of  water  are  lost  in  sinkholes.  He  recommended  trac¬ 
ing  shallow  near-horizontal  ground-water  flows  using  radial  configurations  of 
electrodes  and  reported  that  streaming  potentials  are  a  maximum  in  the  direc¬ 
tion  of  flow. 

13.  Cooper,  Koe*-  er,  and  Franklin  (1982)  reported  that  the  electrochem¬ 
ical  and  telluric  contribution  to  SP  are  relatively  minor  and  that  one  may 
assume  that  values  measured  relative  to  a  local  reference  in  a  similar  geo¬ 
logic  condition  are  due  to  the  electroklnetic  contribution  or  streaming  poten¬ 
tial.  Likewise,  Corwin  and  Hoover  (1978)  have  demonstrated  that,  for  models 
of  similar  geometry,  SP  anomalies  generated  by  electroklnetic  coupling  are 
larger  in  amplitude  than  those  generated  by  other  mechanisms.  Nourbehecht 
(1963)  determined  that  streaming  potential  mechanisms  are  capable  of  producing 
anomalies  of  "several  hundred  millivolts."  These  figures  are  consistent  with 
the  anomalies  commonly  measured  in  the  present  study  (as  great  as  -500  mv) . 
Streaming  potentials  mav  produce  very  large  anomalies  on  the  surface.  Accord¬ 
ing  to  Corry  (1985),  areas  with  large  amplitude  streaming  potentials  are  char¬ 
acterized  by  high  rainfall  and  steep  topography.  An  .xtreme  example  is  the 
-2,693-mv  anomaly  found  on  Agadak  Volcano  in  Alaska  (Corwin  and  Hoover  1978) 
and  the  -1,600-mv  anomaly  associated  with  Kilauea  Volcano  in  Hawaii. 

14.  The  streaming  potentials  produced  by  the  motion  of  ground  water 
under  a  pressure  gradient  through  a  porous  media  result  from  the  preferential 
adsorption  of  ions.  This  process  was  first  recognized  in  1809  in  its  inverse 
form — electro-osmosis,  where  the  application  of  an  electric  potential  across  a 
wet  soil  causes  a  net  flow  of  water  to  the  cathode  (Mitchell  1976).  In  natu¬ 
ral  systems,  streaming  potentials  occur  because  of  the  electric  double  layer 
of  ions  associated  with  the  interface  between  the  mineral  grains  and  the  pore 
fluid  (ground  water) .  Unsatisfied  bonds  at  the  surface  of  mineral  grains 
adsorb  ions  from  solution;  positive  ions  are  attracted  to  the  surface,  nega¬ 
tive  ions  are  repulsed.  The  direction  of  flow,  or  streaming,  is  therefore 
characterized  by  a  net  increase  of  negative  ions  in  the  solution  (Figure  1). 
The  phenomenon  was  explained  and  quantified  by  Helmholtz  in  1879  and  has  been 
further  refined  by  the  work  of  others — especially  soil  physicists  and  colloi¬ 
dal  chemists. 

15.  According  to  the  Helmholtz  equation. 
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where 

V  =  the  measured  stream  potential  (stat  volts) 

stat  volt-sec-ft\ 

p  =  fluid  resistivity  (stat  ohm-ft  or  - — — - - — r — ) 

■'  \  stat  coulomb  / 

^  ^  ^  /stat  coulomb\ 

e  *  fluid  dielectric  constant  ( — — r - rr“rr) 

\stat  volt-ft/ 

g  =  voltage  across  the  Helmholtz  double  layer  (stat  volts) 
y  =  viscosity  of  pore  fluid 

AP  =  pressure  drop  along  flow 

16.  For  a  given  set  of  materials  properties,  the  only  variables  in  the 
above  equation  are  V  and  AP  ;  V  is  proportional  to  AP  and  V/AP  is  a 
constant  (Cooper,  Koester,  and  Franklin  1982).  V/AP  is  known  as  the  electro- 
kinetic  coupling  coefficient,  but  it  is  difficult  to  calculate  because  little 
is  known  about  the  behavior  of  p  ,  g  ,  and  u  in  the  pores  of  rocks  and 
soil  (Corwin  and  Hoover  1978).  For  laminar  flow  through  porous  media,  the 
quantities  V  and  AP  are  related  through  Darcy's  law  to  permeability, 
hydraulic  gradient,  and  flow.  Sand-sized  sediment  produces  the  greatest 
streaming  potential  (i.e.,  large  negative  values).  Streaming  potential  values 
drop  (become  more  positive)  for  media  of  lower  permeabilities  (silts  and 
clays)  and  for  flow  through  coarse  rubble  in  open  fissures  (Ogilvy,  Ayed,  and 
Bogoslovsky  1969)  .  If  turbulent  fluid  flow  takes  place  in  channels  and  frac¬ 
tures,  it  is  thought  that  less  Intense  streaming  potentials  are  caused  than 
would  arise  from  an  equal  quantity  of  water  flowing  through  a  porous  media 
under  laminar  conditions. 

17.  Understanding  the  relationship  between  SP  and  ground-water  flow 
through  fractures  and  channels  is  critical  to  a  successful  application  of  the 
technique  in  areas  of  consolidated  bedrock.  This  is  especially  true  for 
karst,  where  flow  regimes  can  range  from  diffuse  fracture  flows  to  open  solu- 
tional  conduits.  Often,  in  karst,  the  solutlonally  enlarged  channels  are 
wholly  or  partially  filled  with  secondary  material  which  can  consist  of  par¬ 
ticle  sizes  ranging  from  clay  to  gravels.  Karst  channels  can  also  be  either 
air  filled,  water  filled,  or  a  combination  of  air,  water,  and  sediments.  Such 
complexity  requires  careful  consideration  in  interpreting  SP  results  over 
these  areas. 
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18.  The  literature  is  scant  regarding  streaming  potentials  produced  in 
fracture  or  conduit  flow  conditions.  Laboratory  experiments  have  been  con¬ 
ducted  by  Ahmad  (1964)  and  Bogoslovsky  and  Ogllvy  (1972)  to  study  streaming 
potentials  in  channel  flow  models.  Both  studies  affirm  the  proportional  rela¬ 
tionship  between  SP  and  AP  ,  or  hydraulic  head  (the  driving  force  behind  the 
fluid  movement).  Otherwise,  they  found  that  SP  fields  in  fissured  media  dif¬ 
fer  from  those  in  unconsolidated  porous  media.  For  Instance,  at  similar  pres¬ 
sures  (AP)  the  SP  values  produced  in  fissures  2  to  6  mm  in  width  are  several 
orders  of  magnitude  smaller  than  those  produced  in  media-grained  sand 
(Bogoslovsky  and  Ogilvy  1972).  SP  values  were  found  to  decrease  with  the 
opening  of  fissures,  l.e.,  larger  fissures  (6  mm)  produce  SP  values  approach¬ 
ing  zero  or  positive  values  while  smaller  fissures  (2  mm)  produce  SP  values 
that  are  more  negative.  However,  the  presence  of  sand  fill  in  the  fissures 
considerably  increased  the  SP  over  values  obtained  for  empty  fissures.  A 
filling  of  40  percent  sand  produced  the  maximum  SP;  whereas,  introducing  a 
clay  fraction  to  the  fill  resulted  in  a  decrease  in  the  SP.  Ahmad  (1964)  con¬ 
cluded  that  the  geometry  of  flow  was  relatively  unimportant  in  its  effect  on 
SP  in  fissured  flow  and  that  SP  decreases  as  salinity  of  the  electrolyte  (pore 
fluid)  Increases. 

19.  Ogllvy,  Ayed,  and  Bogoslovsky  (1969)  reported  that  clay-filled  fis¬ 
sures  with  low  permeabilities  are  characterized  by  low  SP  values.  As  support¬ 
ing  evidence  for  this,  they  cited  the  relatively  high  positive  values  of  SP 
associated  with  clay  layers  in  borehole  geophysical  logs.  They  also  claimed 
that  salinization  also  leads  to  a  decrease  in  SP  values.  Warriner  and  Taylor 
(1986)  conducted  laboratory  experiments  in  an  attempt  to  model  streaming 
potential  mechanisms  and  parameters.  They  also  found  that  the  presence  of 
clay  in  or  near  the  flow  path  degrades  SP  magnitudes,  and  also  mentioned  the 
positive  SP  values  for  shale  layers  in  well  logs  as  confirmation.  In  their 
experiments,  all  measured  streaming  potentials  were  negative. 

20.  The  dependency  between  SP  values  and  AP  has  been  further  demon¬ 
strated  at  various  Corps  of  Engineers  reservoirs.  Perturbations  to  the  water 
flow  regimes  were  created  by  changing  the  level  of  the  reservoir  pool.  These 
changes  in  hydraulic  head  correlated  to  SP  fluctuations  at  electrodes  located 
over  seepage  zones,  indicating  that  flow  rate  affects  SP  (Warriner  and  Taylor 
1982). 
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21.  Some  insight  into  depth  resolution  and  the  ability  of  SP  to  detect 
water-filled  conduits  was  gained  in  an  SP  survey  over  Manatee  Springs  in 
Florida  (Cooper  1983).  Manatee  Springs  is  a  large  water-filled  solution  con¬ 
duit  and  was  detected  by  SP  at  a  depth  of  90  ft.*  Florida's  karst  often  has 
high  diffuse  permeability  and  may  also  contain  appreciable  amounts  of  sand; 
this  may  explain  the  significant  streaming  potential  anomalies  encountered  at 
this  test  site.  Further  evidence  of  SP  depth  capability  was  demonstrated  by  a 
Corps  of  Engineers  survey  at  Waterbury  Dam  in  Vermont,  where  there  was  a  good 
correlation  of  negative  anomalies  with  man-made  outlet  works  conduits  at  a 
depth  of  120  ft  below  the  embankment  on  which  the  surveys  were  conducted. 

22.  It  is  important  to  refine  the  SP  technique  and  our  interpretation 
of  its  results  to  distinguish  between  different  types  of  ground-water  flow 
regimes.  Most  basic  would  be  dif ferentation  between  filtration  flow  (seepage) 
through  soil  and  unconsolidated  sediments  and  concentrated  flow  in  fractures 
or  channels.  Within  karst  terrain  this  amounts  to  detecting  shallow  recharge 
ground  water  percolating  through  soil  and  regolith,  water  confined  to  solu- 
tionally  developed  channels  in  the  eplkarst  and  moving  laterally  or  verti¬ 
cally,  and  deeper  flow  through  solution  conduits  in  either  vadose  or  phreatic 
conditions.  This  variety  of  conditions  entails  a  multitude  of  variables  which 
could  affect  SP  generation  and  interpretation,  such  as: 

a.  Local  relief,  which  could  influence  hydraulic  gradients. 

b.  Lithology  (insolubles  would  affect  permeability  of  fill). 

£.  Variations  in  water  and  soil  chemistry. 

d.  Bedrock  structure  and  deformation. 

e.  Degree  of  karstlf Ication,  extent  of  subsurface  drainage. 

f.  Depth  to  "water  table"  or  major  conduits. 

£.  Climatic  and  seasonal  parameters. 


*  A  table  of  factors  for  converting  non-SI  to  SI  (metric)  units  of  measure¬ 
ment  is  presented  on  page  3. 
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PART  II:  DESCRIPTION  OF  STUDY  AREA 


Sites 


23.  In  this  research,  SP  was  studied  by  means  of  observations  and  data 
from  a  number  of  sites  in  Rockbridge  County,  VA — within  the  Valley  and  Ridge 
Province.  All  of  the  sites  are  underlain  by  Ordovician  limestones  and  dolo¬ 
mites.  Karst  topography  is  well  developed  locally;  numerous  sinkholes,  caves, 
karren,  and  sinking  streams  are  present.  Elevation  is  about  360  m  above  sea 
level. 

24.  Nine  sites  were  selected  for  study.  Six  sites  are  sinkholes 

(one  contains  a  cave  entrance),  two  sites  are  sinking  streams,  and  one  site  is 
at  an  intermittent  stream  bed.  Figure  2  shows  the  location  of  Rockbridge 
County  and  several  of  the  study  sites.  The  sites  were  chosen  to  represent  a 
range  of  karst  features  (collapse  sinks,  subsidence  sinks,  sinkholes  with 
known  cave  passages  under  them,  sinkholes  above  springs,  fault-related  sink¬ 
holes  and  sinking  streams,  and  a  failed  sinkhole  pond).  All  are  examples  of 
ground-water  recharge  points,  and,  as  such,  should  exhibit  high  streaming 
potential  values. 

Harris  and  Hunter  sinkholes 

25.  Two  of  the  sites  (Harris  and  Hunter  sinkholes)  were  monitored  in 
1985;  for  all  the  other  sites,  study  began  in  July  1986.  Unfortunately,  the 
commencement  of  this  research  in  1986  coincided  with  a  prolonged  and  severe 
drought  that  caused  record  low  water  levels  throughout  the  area.  The  natural 
oscillations  in  water  levels  and  flow  volumes  that  were  anticipated  in  the 
design  of  this  research  did  not  occur  until  the  spring  of  1987.  The  Harris 
and  Hunter  sinkholes  are  located  in  open  pasture  land  on  Kiger  Hill  and  are 
developed  in  the  Beekmantown  Formation  on  an  anticlinal  fold.  The  Harris 
sinkhole  collapsed  during  heavy  rain  in  August  1984  and  now  displays  a  blind 
shaft  at  its  bottom.  The  Hunter  sinkhole  is  a  more  shallow  subsidence  feature 
and  is  located  180  m  up-gradient  from  the  Harris  site. 

Moore  sinkholes 

26.  The  Moore  sinkhole  site  is  also  in  open  pasture  land  about  1  mile 
west  of  Lexington.  It  is  underlain  by  the  Lincolnshire  limestone  which  dis¬ 
plays  numerous  large  sinkholes  in  the  immediate  area.  The  Moore  sinkhole  "A" 
was  developed  into  a  farm  pond  in  the  1970 's  and  apparently  held  water 
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I  (surface  runoff)  successfully  for  several  years  until  leakage  failure  occurred 

,  around  1980.  It  has  been  dry  since  then,  although  there  is  currently  another 

pond  holding  water  90  m  away  (which  is  not  a  sinkhole).  The  Moore  sinkhole  A 
is  aligned  with  other  adjacent  sinkholes  in  a  northeasterly  bearing  (along 
strike) .  A  small  transverse  fault  (striking  E-W)  is  exposed  one-half  mile 
east  of  the  site.  Two  other  sinkholes  (Moore  "B"  and  "C")  are  located  down- 
gradient  from  and  close  to  sinkhole  A.  These  three  sinkholes  form  a  related 
complex  and  were  monitored  for  SP,  electrical  resistivity  (ER) ,  and  electro¬ 
magnetic  conductivity  (EM)  data  beginning  in  December  1986  to  determine  the 
nature  of  their  relationship  and  to  look  for  any  hydrologic  connections 
implied  by  the  data. 

McHenry  sinkholes  and  cave 

27.  The  McHenry  site  contains  two  sinkholes  (A  and  B)  and  a  cave  named 
Whose  Hole.  This  area  is  located  at  the  base  of  Little  Camp  Mountain  (local 
relief  is  450  m)  in  the  Lincolnshire  limestone.  The  structure  here  is  com¬ 
plex,  consisting  of  a  series  of  small,  compressed  folds  and  minor  faulting. 
Whose  Hole  Cave  contains  a  permanent  stream  which  flows  northeast,  along 
strike,  for  0.4  km  to  discharge  at  a  spring  immediately  below  sinkhole  B. 
Sinkhole  B  is  developed  in  the  trough  of  a  small  syncline  and  contains  a  rela¬ 
tively  recently  collapsed  blind  shaft  (although  no  connection  can  be  made  here 
into  the  cave  passage  that  is  presumed  to  underlie  it) .  A  fluorescent  dye 
trace  established  the  hydrologic  connection  between  the  cave  stream  and  the 
spring  at  sinkhole  B.  Sinkhole  A  at  the  McHenry  site  is  located  between  the 
cave  and  sinkhole  B.  It  is  a  shallow  but  well  defined  feature  and  contains 
good  exposures  of  bedrock  that  suggest  that  a  fault  cuts  through  this  sink¬ 
hole.  An  intermittent  stream  course  is  adjacent  to  the  sinkhole  entrance  of 
Whose  Hole  Cave  and  trends  easterly  to  join  the  stream  that  flows  from  the 
spring  discharge  of  the  cave.  A  long-term  electrode  monitoring  station  has 
been  Installed  across  this  normally  dry  streambed  at  a  point  270  m  east  of  the 
cave  entrance.  All  of  the  McHenry  sites  are  in  pasture  land  except  for  the 
cave  site  which  is  wooded. 

Spring  Branch  sinking  stream 

28.  One  of  the  sinking  stream  sites  is  located  on  Spring  Branch  about 
1.6  km  east  of  the  McHenry  site.  The  Spring  Branch  site  is  also  developed  in 
the  Lincolnshire  limestone.  This  is  technically  a  loosing  stream  in  that  nor¬ 
mally  only  a  portion  of  the  flow  in  Spring  Branch  is  lost  underground  here 


3 


through  a  series  of  four  or  five  enlarged  joints  in  the  bed  and  bank  of  the 
stream.  However,  between  July  1986  and  January  1987,  the  entire  flow  of 
Spring  Branch  was  lost  underground  because  of  the  low  flow  levels  resulting 
from  drought.  A  thrust  fault  runs  parallel  to  Spring  Branch  and  has  been  sug¬ 
gested  as  a  plausible  cause  of  the  watershed  piracy  that  has  been  observed  in 
the  area  (Johnson  1980  and  MacMillan  1980) .  It  is  likely  that  the  Spring 
Branch  site  is  related  to  this  fault  and  that  the  water  lost  from  Spring 
Branch  is  diverted  to  another  watershed  (South  Buffalo  Creek). 

Ford  site's  sinking  stream 

29.  The  Ford  site  is  another  example  of  a  sinking  stream;  overflow  from 
a  spring  is  lost  into  solutlonal  fissures  in  a  rocky  streambed  1  km  west  of 
Lexington.  This  site  is  underlain  by  Lincolnshire  limestone  as  well.  The 
extreme  surface  rockiness  and  numerous  sinkholes  nearby  suggest  a  high  degree 
of  solutlonal  weathering  has  taken  place  here.  During  very  heavy  precipita¬ 
tion  some  of  the  water  continues  beyond  the  normal  point  of  loss  for  about 
45  m  to  the  entrance  of  a  small  cave  passage  where  it  is  finally  lost  to 
underground  channels. 


Research  Methodology 

30.  An  array  of  electrodes  was  Installed  at  each  site  and  monitored  for 
SP  over  a  period  of  months.  A  double  ring  concentric  configuration  was  used 
for  sinkholes  and  a  rectilinear  grid  pattern  was  used  at  sinking  stream  sites 
and  as  a  tracking  array  between  sinkholes.  Electrodes  were  numbered  with  ref¬ 
erence  to  the  compass  so  that  spatial  patterns  of  anomalies  could  be  related 
to  other  directional  features  (such  as  geologic  structure,  stream  flow  direc¬ 
tions,  etc.).  A  reference  electrode  was  located  at  least  25  m  distant  from 
the  test  electrodes  in  a  direction  up-surface  elevation  gradient  from  the 
array.  Some  sites  received  two  or  more  reference  electrodes  (for  Instance,  in 
the  center  of  a  sinkhole)  as  a  check  on  the  results. 

31.  SP  values  were  read  for  each  site  at  approximately  weekly  inter¬ 
vals,  with  special  attention  following  precipitation  events.  Although  SP  val¬ 
ues  for  electrodes  varied  considerably  over  time,  the  relationship  between 
electrodes  remained  constant  (i.e.,  the  strongly  negative  anomalies  were 
always  in  the  same  location) .  The  locations  of  consistently  anomalous  elec¬ 
trodes  were  presumed  to  represent  areas  of  concentrated  ground-water  flow 
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(maximum  streaming  potential).  SP  readings  were  taken  with  a  high  impedance 
digital  voltmeter  by  connecting  wire  leads  between  the  reference  electrode 
(always  common,  or  ground,  for  polarity  convention)  and  each  test  electrode. 

SP  values  ranged  between  +400  and  -500  mv  and  were  steady  while  reading  (fluc¬ 
tuating  only  one  or  two  mv) . 

32.  Other  types  of  data  were  collected  to  confirm  interpretation  of  SP 
anomalies  as  areas  of  increased  ground-water  flow.  Electrical  resistivity 
(both  Wenner  and  Lee  arrays)  and  electromagnetic  conductivity  surveys  were 
conducted  at  the  Moore  and  cave  sites  to  provide  an  alternate  geophysical 
interpretation  of  parameters  such  as  anisotropy  and  permeability.  Precipita¬ 
tion  (from  rain  gages)  and  soil  temperature  data  (by  thermometer  and  thermo¬ 
couple  probes)  were  recorded  for  each  site.  Fluctuations  in  stream  levels 
were  determined  by  constructing  weirs  above  sinking  stream  sites  and  from  a 
float/electrical  resistance  device  in  the  cave  stream.  In  addition,  visual 
observation  of  surface  runoff  patterns  along  with  geological  study  of  nearby 
outcrops  provided  clues  for  interpreting  subsurface  conditions  and  related  SP 
effects . 

33.  The  electrodes  used  were  copper-clad  steel  rods,  2  cm  in  diameter 
and  60  cm  in  length,  and  were  driven  into  the  ground  about  50  cm.  This  type 
of  electrode  was  successfully  used  to  detect  SP  anomalies  over  zones  of  leak¬ 
age  In  dams  by  the  Corps  of  Engineers  (results  were  confirmed  by  both  drill¬ 
ing  and  other  geophysical  methods).  This  Is  a  departure  from  the  traditional 
technique  of  using  nonpolarizing  (NP)  SP  electrodes.  NP  electrodes  are  com¬ 
monly  employed  in  mineral  and  geothermal  surveys  and  in  corrosion  control 
applications.  NP  electrodes  consist  of  a  copper  rod  immersed  in  a  saturated 
solution  of  copper  sulphate  with  a  porous  ceramic  tip  that  allows  ionic  dif¬ 
fusion  and  electrical  contact  between  the  soil  and  electrolyte.  The  purpose 
of  this  design  is  to  eliminate  the  effects  of  local  polarization  (concentra¬ 
tion  cells)  that  may  occur  on  metallic  electrode  surfaces. 

34.  Although  the  use  of  metal  electrodes  does  introduce  some  local 
polarization  effects  into  SP  data.  It  appears  that  they  are  effective  in  long¬ 
term  monitoring  of  streaming  potential  and  that  polarization  effects  are  mini¬ 
mal  and  tend  to  be  averaged  out.  The  NP  electrodes  would  be  preferable  for 
Instant  readings  of  one-time  surveys  over  large  areas  and  for  applications 
where  anomalies  are  small  and  accuracy  requirements  are  great.  For  use  in 
small  areas  of  uniform  geologic  setting,  metal  electrodes  are  apparently 
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adequate  for  long-term  monitoring  of  large-scale  anomalies  and  systematic 
changes  in  SP.  A  process  of  equilibrium  occurs  for  metal  electrodes  after 
being  placed  in  the  ground;  the  time  required  for  this  may  vary  from  1  day  to 
1  week.  SP  values  become  relatively  stable  when  this  equilibrium  of  polari¬ 
zation  effects  takes  place.  The  polarization  equilibrium  process  for  several 
test  site  electrodes  is  recorded  in  Figure  3.  Adequate  time  should  be  allowed 
before  reading  SP  values  after  initial  electrode  placement.  Electrodes  are 
then  left  in  place  for  as  long  as  needed.  Polarization  effects  are  further 
minimized  because  the  high  Internal  resistance  of  the  metres  that  are  used 
allows  practically  no  current  flow  to  occur.  The  use  of  a  bimetallic  elec¬ 
trode  (copper-clad  steel)  may  be  advantageous  in  that  the  exposed  steel  tip 
acts  as  a  sacrificial  anode,  thereby  reducing  corrosion  and  polarization 
effects  on  the  copper  sheath.  The  advantages  of  metal  electrodes  over  the  NP 
electrodes  are  their  low  cost,  ease  of  placement  and  use,  and  durability.  NP 
electrodes  were  used  to  verify  the  results  obtained  from  metal  electrodes  at 
some  of  the  sites  in  this  research.  The  actual  SP  values  were  different,  but 
in  most  cases  there  was  general  agreement  of  relative  values  between  electrode 
locations  (lowest  values  were  in  the  same  places,  etc.). 

35.  In  addition  to  the  manual  collection  of  data  with  voltmeters,  an 
automated  data  collection  system  was  devised  and  Implemented.  This  system 
consisted  of  a  portable  battery  powered  computer,  an  analog/digital  signal 
converter  Interface  and  temperature/soil  moisture  sensors.  The  entire  system 
was  housed  in  a  weatherproof  case  and  could  be  left  in  the  field  and  pro¬ 
grammed  to  gather  data  at  any  desired  frequency.  The  metal  electrodes  were 
wired  directly  to  the  signal  Inputs  of  the  Interface.  The  automated  data  col¬ 
lection  (ADC)  system  was  used  for  two  tasks — to  monitor  the  effects  of  envi¬ 
ronmental  variables  on  SP  over  time  and  to  monitor  SP  changes  induced  by  the 
artificial  recharge  (dumping)  of  water  into  a  sinkhole.  The  long-term  moni¬ 
toring  of  electrode  behavior  was  conducted  at  the  intermittent  streambed  site 
east  of  the  cave  on  the  McHenry  site.  The  water  dumping  experiment  took  place 
at  the  Moore  sinkhole  A  near  Lexington.  An  example  of  the  ADC  record  from  the 
long-term  monitoring  is  shown  in  Table  1,  and  from  the  water  dumping  test  in 
Table  2. 

36.  Two  water  dump  tests  were  conducted — 4  and  18  October  1986.  For 
both  tests  15,900  SL  of  water  were  discharged  into  the  center  of  the  sinkhole 
from  four  trucks  over  a  period  of  several  minutes.  SP  measurements  were  taken 


16 


before,  during,  and  after  the  dumping,  both  by  voltmeter  and  by  the  ADC  sys¬ 
tem.  The  purpose  of  this  test  was  to  look  for  SP  response  to  the  flow  of  a 
known  amount  of  water  into  the  sinkhole  system,  with  the  hope  that  SP  changes 
could  be  related  to  flow  rate  in  a  soll-regolith  recharge  setting. 

37.  Other  experiments  aimed  at  understanding  SP  changes  relative  to 
ground-water  flow  were  conducted  at  Whose  Hole  Cave  on  the  McHenry  site.  The 
cave  was  accurately  surveyed,  and  electrodes  were  placed  on  the  surface  at 
locations  that  could  be  related  to  underlying  cave  features.  A  weir  was  con¬ 
structed  across  the  cave  stream  and  a  flow  measuring  device  was  installed  and 
wired  to  the  surface  so  that  stream  level  changes  could  be  conveniently  moni¬ 
tored.  This  device  consisted  of  a  float  attached  to  a  potentiometer;  changes 
in  float  height  could  then  be  detected  as  a  linear  change  of  resistance  and 
calibrated  in  relation  to  water  levels  at  the  weir. 

38.  An  array  of  electrodes  was  installed  on  the  surface  In  rows  perpen¬ 
dicular  to  the  cave  stream  passage  at  a  point  directly  above  the  weir  loca¬ 
tion.  The  vertical  distance  from  the  center  of  this  array  to  the  cave  stream 
was  calculated  to  be  25  m.  In  addition,  a  ring  of  electrodes  was  placed 
around  the  sinkhole  entrance  to  the  cave  so  that  SP  data  could  be  related  to 
underlying  cave  features  at  depths  ranging  from  10  to  20  m.  An  initial  ring 
of  12  electrodes  was  used  from  September  to  December,  1986,  after  which  the 
standard  double  ring  of  32  electrodes  was  installed  in  the  sinkhole.  In 
December  1986  an  EM  conductivity  survey  was  conducted  at  this  sinkhole  using 
Intercoil  spacings  of  10  m. 

39.  For  all  the  sites,  the  SP  data  show  considerable  variability  as 
seen  in  the  standard  deviations  of  individual  electrode  readings.  Although 
the  relative  values  between  electrodes  remain  fairly  constant,  it  is  difficult 
to  interpret  meaningful  trends  from  masses  of  variable  numbers.  Others  have 
attempted  to  gain  better  resolution  in  their  data  by  determining  and  removing 
a  mean  regional  SP  value  (Cooper  1983) .  The  mean  regional  SP  is  determined  by 
averaging  all  the  SP  values  for  a  given  reading  and  subtracting  this  value 
from  each  individual  value  for  that  date. 

40.  The  present  research  suggests  that  the  time  varying  influence  of 
temperature  and  moisture  introduces  systematic  shifts  in  the  SP  data.  This 
influence  can  mask  or  exaggerate  the  streaming  potential  values.  An  effort 
was  made  to  filter  the  data  by  using  ranks  to  preserve  the  ability  of  the  data 
to  document  negative  anomalies  and  variability  while  removing  the  seemingly 
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erratic  swings  In  the  raw  data.  Negative  readings  are  shown  as  low  rankings 
and  positive  readings  are  shown  as  high  rankings  (Erchul  1987).  The  major 
advantage  Is  the  removal  of  systematic  changes  In  the  data.  The  major  disad¬ 
vantage  Is  the  Insensitivity  to  factors  which  exaggerate  scale.  A  more  nega¬ 
tive  reading  due  to  Increased  ground-water  flow  is  not  visible  in  the  ranked 
data  until  one  reading  passes  another.  A  comparison  of  the  spatial  patterns 
of  SP  anomalies  calculated  by  both  averaging  and  ranking  shows  that  no  signif¬ 
icant  difference  is  obtained.  This  suggests  that  the  use  of  averaging  may  be 
adequate  for  plotting  SP  values  to  infer  flow  paths  and  that  systematic  data 
changes  of  whatever  source  are  cancelled  by  averaging.  Further  refinements  in 
quantifying  SP  data  for  use  In  flow  rate  studies,  however,  will  probably 
necessitate  development  of  additional  statistical  treatment. 
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PART  III:  RESULTS  AND  DISCUSSION 


Evaluating  and  Comparing  SP  Results 

41.  This  research  set  out  to  examine  several  basic  questions  by  evalu¬ 
ating  and  comparing  SP  results  from  various  sites  having  different  character¬ 
istics.  Of  primary  concern  are  the  following  issues: 

a.  How  effective  is  SP  in  detecting  flow  paths  in  karst  features? 

b.  Are  detected  anomalies  characterized  by  percolation  and  seepage 
flow  or  by  fracture-conduit  flow? 

£.  Can  SP  results  be  related  to  flow  rates? 

d.  What  is  the  relationship  between  SP  variations  and  environ¬ 
mental  variables? 

42.  At  every  site  significant  SP  anomalies  have  been  measured.  These 
anomalies  differ  by  as  much  as  an  order  of  magnitude  from  the  SP  values  of  the 
surrounding  electrodes,  and  their  locations  remain  fixed  within  the  electrode 
array  space.  The  magnitudes  vary  with  time,  but  the  relative  values  between 
electrodes  remain  fairly  consistent.  At  the  Harris  and  Hunter  sinkholes, 
anomalies  were  linked  and  used  to  infer  flow  paths  within  and  between  the  two 
sinkholes  over  a  distance  of  180  m  (Erchul  1987).  The  plotted  SP  values  were 
calculated  by  averaging  all  readings  for  each  date  and  also  by  ranking  the 
data  (lowest  reading  electrode  for  a  given  date  equal  to  one,  etc.).  For  pur¬ 
poses  of  plotting  the  values  to  infer  flow  paths,  the  spatial  relationships 
produced  by  these  two  methods  of  calculation  do  not  differ  significantly 
(i.e.,  the  anomalies  Ho  not  change  relative  to  each  other,  and  connecting  them 
creates  the  same  flow  paths).  At  the  Harris  and  Hunter  site,  the  flow  paths 
are  aligned  with  the  geologic  strike  of  the  underlying  bedrock  (N  50®  E)  and 
with  joint  or  fracture  orientations  perpendicular  to  the  strike.  An  electri¬ 
cal  resistivity  survey  showed  anisotropy  in  the  form  of  an  ellipse  oriented  in 
the  strike  direction,  also.  In  addition,  nearby  caves  (there  are  five  within 
1.5  km  of  the  site)  are  predominately  strike  controlled  in  their  passage 
development  (Figure  4).  It  is  reasonable  to  assume  that  the  flow  paths 
inferred  by  the  SP  data  at  the  Harris  and  Hunter  site  represent  strike  and 
offset  joint  control  that  has  created  areas  of  preferential  water  movement. 
Figure  5  shows  electrode  arrays  and  inferred  flow  paths  at  the  Harris  and 
Hunter  sinkholes  (Erchul  1987) . 
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Spring  Branch  site  results 

43.  The  Spring  Branch  site  was  selected  because  it  is  a  sinking  stream 

an  ’  appears  to  be  related  to  a  fault  that  may  divert  ground  water  out  of  the 
watershed.  A  grid  array  of  electrodes  spaced  on  7.5-m  centers  was  placed  on 

the  hillside  above  the  point  of  stream  loss  (swallet).  A  weir  was  constructed 

above  the  swallet  to  measure  fluctuations  in  flow.  An  experiment  was  con¬ 
ducted  with  two  reference  electrodes  at  this  site.  The  results  are  shown  in 
Figure  6  and  indicate  excellent  agreement  between  SP  values  for  all  elec¬ 
trodes.  Data  were  collected  weekly  for  about  3  months  and  used  to  calculate 

average  SP  values.  Figure  7  shows  the  Spring  Branch  site  with  the  electrode 

array  and  an  SP  (mv)  contour  map.  Anomalous  areas  are  identifiable  as  concen¬ 
trations  of  contours  (areas  of  maximum  SP  gradient).  The  bedrock  here  strikes 
northeast  and  dips  15°  SE.  The  stream  sinks  into  a  solutlonally  enlarged 
joint  striking  SE  (direction  of  dip)  which  coincides  with  the  orientation  of 
the  zone  of  maximum  SP  gradient.  The  contours  indicate  preferential  water 
movement  in  this  southeasterly  direction  and  also  to  the  east.  From  examina¬ 
tion  of  outcrops  exposed  along  the  streambed,  a  set  of  joints  has  been  mea¬ 
sured  striking.  N  85°  E.  The  two  linear  trends  displayed  on  the  SP  contour  map 
and  marked  B-B  and  C-C  are  interpreted  as  ground-water  flow  paths  based  on 
observation  of  outcrop  features  (dip  and  joint  strikes  being  the  dominant 
structural  controls).  No  physical  evidence  of  faulting  is  displayed  at  the 
site,  but  the  geologic  map  (Spencer  1968)  indicates  a  low  angle  thrust  fault 
in  the  vicinity,  dipping  east.  MacMillan  (1980)  and  Johnson  (1980)  have  iden¬ 
tified  this  fault  and  this  portion  of  Spring  Branch  as  the  most  likely  source 
of  watershed  loss  into  the  adjacent  watershed  to  the  east.  The  SP  data  appear 
to  confirm  this  interpretation  by  indicating  a  major  anomalous  zone  below  the 
swallet.  This  could  be  the  point  where  the  fault  captures  the  water  lost  from 
Spring  Branch. 

44.  Since  the  water  from  Spring  Branch  is  lost  into  solutlonal  chan¬ 
nels,  this  site  should  provide  some  indication  of  SP  ability  to  detect  chan¬ 
nelized  flow.  Averages  were  calculated  for  the  western  half  of  the  electrode 
array  (upstream  of  the  swallet)  and  for  the  eastern  half  (downstream  of  the 
swallet).  The  results  are  shown  in  Figure  8,  along  with  a  record  of  stream 
flow  and  precipitation.  The  SP  values  for  the  eastern  half  of  the  array  are 
more  negative  than  the  western  half,  indicating  that  those  electrodes  down- 
gradient  of  the  swallet  may  be  measuring  increased  streaming  potential  because 
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of  channelized  flow  below  them.  The  entire  array  is  located  on  a  uniformly 
north  sloping  hillside  with  no  obvious  differences  between  the  two  halves.  It 
can  also  be  seen  from  Figure  8  that  SP  becomes  more  negative  in  response  to 
either  precipitation  (soil  moisture)  and/or  streamflow.  Figure  9  depicts  evi¬ 
dence  for  SP  generation  because  of  the  nonchannellzed  flow  at  Spring  Branch 
for  those  electrodes  nearest  the  streambed.  The  SP  data  for  the  two  dates 
shown  indicate  a  negative  shift  of  about  200  mv  at  electrodes  Nos.  31,  4a,  and 
5a  during  a  period  of  increased  streamflow.  These  electrodes  are  only  about 
5  m  from  the  bed  of  the  stream.  On  10  October  the  streambed  below  the  swallet 
was  dry  (all  the  flow  was  going  into  the  swallet)  but  on  2  November  the  added 
flow  caused  the  swallet  to  overflow  and  this  additional  water  was  absorbed 
into  the  sediments  of  the  streambed  for  a  distance  of  25  m  below  the  swallet. 
This  seepage  flow  through  the  streambed  sediments  probably  caused  an  increased 
streaming  potential  that  is  represented  by  the  negative  shift  in  those 
streamslde  electrodes  as  depicted  in  Figure  9.  A  similar  negative  SP  response 
was  measured  at  an  electrode  of  the  long-term  monitoring  site  whenever  the 
normally  dry  streambed  would  carry  water  following  heavy  rains. 

Ford  site  results 

45.  The  Ford  site  is  a  sinking  stream  also,  although  no  faulting  is 
known  to  occur  here.  Overflow  from  a  spring  development  for  livestock,  along 
with  intermittent  surface  water  runoff,  is  lost  into  solutional  openings  at 
the  head  of  a  small  steep-sided  valley.  Flood  flow  produces  more  water  than 
can  be  taken  by  the  swallet  and  the  excess  is  carried  25  m  beyond  to  the 
entrance  of  a  small  cave  (a  single  conduit  7  m  in  length).  Figure  10  shows  a 
map  of  the  Ford  site  with  electrode  array  and  SP  values  plotted  as  millivolt 
contours.  Data  for  this  site  does  not  seem  to  depict  a  linear  flow  path.  The 
greatest  anomaly  is  on  a  line  connecting  the  swallet  and  overflow  cave,  and 
this  line  coincides  with  the  direction  of  strike  (northeast).  Streamflow 
(from  a  weir  above  the  swallet),  soil  temperature,  and  average  SP  for  the 
array  are  recorded,  for  the  period  June  -  November,  along  with  SP  for  elec¬ 
trode  No.  9  (the  greatest  anomaly)  in  Figure  11.  There  is  no  obvious  cor¬ 
relation  here  between  streamflow  and  SP.  Soil  temperature  has  a  direct  effect 
on  array  average  SP  at  times  (the  month  of  August,  especially),  but  at  other 
times  this  relationship  seems  less  direct,  probably  the  result  of  precipita¬ 
tion  events  interacting  with  temperature  change.  The  response  of  the  anoma¬ 
lous  electrode  No.  9  differs  from  the  response  of  the  array  average. 
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indicating  that  No.  9  reacts  differently  to  environmental  changes  or  to  other 
unknown  factor (s)  that  may  Influence  this  anomalous  zone  but  not  the  surround¬ 
ing  electrodes.  It  may  be  that  the  anomaly  at  electrode  No.  9  represents  a 
site  where  ground  water  is  descending  vertically  rather  than  flowing  later¬ 
ally.  Vertical  shafts  are  not  unusual  in  karst  aquifers  and  are  a  principal 
element  in  the  epikarst  model  (subcutaneous  zone)  of  Williams  (1983)  and  White 
(1977).  It  is  reasonable  to  assume  the  existence  of  vertical  ground-water 
movement  at  the  Ford  site  because: 

a.  The  base  level  (Maury  River)  is  40  m  lower  than  the  Ford  site. 

b.  There  are  no  known  major  springs  or  resurgences  between  the 
swallet  and  the  Maury  River  (a  distance  of  2.4  km). 

£.  The  lithology  of  the  site,  and  terrain  between  the  site  and  the 
Maury  River,  is  soluble  carbonates  with  no  perching  beds  or 
aqulcludes . 

46.  An  examination  of  the  SP  data  collected  at  both  of  the  sinking 
stream  sites  (Tables  3  and  4)  shows  an  interesting  seasonal  trend  which  is  not 
evident  at  any  of  the  other  study  sites.  At  both  the  Spring  Branch  and  Ford 
sites  the  SP  values  trend  from  largely  negative  in  the  summer  and  fall  (July 
to  November)  to  almost  entirely  positive  for  the  period  December  through 
April.  The  period  July  to  November  corresponds  to  the  severe  drought  which 
began  to  cease  with  the  return  of  significant  rainfall  in  December.  The  only 
other  sice  to  show  a  pronounced  reaction  to  the  drought  was  the  sinkhole 
entrance  to  Whose  Hole  Cave,  but  here  the  SP  response  was  exactly  opposite 
(positive  during  the  drought  and  negative  in  response  to  returning  normal  rain 
and  ground-water  flow) .  The  fact  that  the  sinking  stream  sites  show  this 
apparently  unique  reaction  may  indicate  that  conduit  flow  is  playing  a  role  in 
SP  values  here,  but  we  do  not  have  an  explanation  for  this  phenomenon  at  this 
time.  It  is  difficult  to  say  why  the  transition  from  a  severe  drought  to  nor¬ 
mal  "wetter"  weather  would  cause  such  wholesale  systematic  shifts  in  SP  from 
negative  to  positive  values. 

McHenry  site  results 

47.  The  McHenry  site  contains  two  sinkholes,  a  cave  with  a  stream  in 
it,  a  spring  (resurgence  of  the  cave  stream),  and  a  typical  dry  karst  (inter¬ 
mittent)  surface  stream  course.  All  of  these  features  are  wichln  0.4  km  of 
each  other  in  an  area  of  complex  geologic  structure  (Figure  12).  The  McHenry 
sinkholes  A  and  B  were  monitored  for  SP  changes  from  July  to  November 
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(Tables  5  and  6).  Figure  12  shows  the  spatial  relationship  between  the  two 
sinkholes  and  the  spring.  Sinkhole  A  is  located  about  15  m  higher  than  sink¬ 
hole  B,  In  a  direction  generally  up-dip  (to  the  southeast).  The  spring 
emerges  from  the  base  of  a  steep  slope  about  10  m  lower  than  sinkhole  B.  The 
flow  of  ground  water  can  be  assumed  to  be  from  A  to  B  to  the  spring.  SP  data 
for  sinkhole  A  is  shown  as  millivolt  contours  in  Figure  13.  The  strongly 
anomalous  zone  in  the  southeast  quadrant  corresponds  to  a  fault  that  is 
exposed  in  bedrock  ledges  outcropping  here.  The  linear  trend  of  these  con¬ 
toured  anomalies  is  oriented  to  the  northwest  in  the  direction  of  sinkhole  B. 
In  Figure  14,  the  fault-related  anomalous  zone  in  sinkhole  A  is  shown  to  have 
an  equal  effect  on  both  the  inner  and  outer  rings  of  electrodes  implying  that 
at  least  some  flow  is  lateral  (towards  sinkhole  B  and  the  spring).  The  milli¬ 
volt  contour  map  for  sinkhole  B  (Figure  15)  exhibits  two  main  anomalous  zones. 
The  anomaly  In  the  southeast  quadrant  probably  represents  the  flow  path  of 
water  coming  from  the  direction  of  sinkhole  A.  The  largest  anomalous  zone,  in 
the  northeast  quadrant,  represents  the  outflow  of  water  from  sinkhole  B  toward 
the  spring.  This  assumption  Is  supported  by  observation  of  outcrops;  the 
spring  emerges  from  the  trough  of  a  small  synclinal  fold.  Ground  water  in  the 
vicinity  of  sinkhole  B  would  therefore  flow  down-dip  on  either  limb  of  the 
syncline  and  would  have  to  exit  along  the  axis  of  the  synclinal  trough.  The 
SP  results  from  McHenry  sinkholes  A  and  B  exhibit  especially  strong  anomalous 
zones  that  demonstrate  a  linear  pattern  of  connections.  The  fortuitous  expo¬ 
sure  of  outcrops  in  the  vicinity  helps  to  confirm  these  interpretations  of 
flow  paths. 

48.  During  the  course  of  monitoring  these  sinkholes,  the  phenomena  of 
short-term  anomalies  were  observed  at  both  sinkholes.  In  both  cases,  the 
events  occurred  between  mid-August  and  mid-October.  Figure  16  shows  an  SP 
Increase  of  500  mv  for  two  electrodes  at  sinkhole  A  (these  two  electrodes  were 
unique  in  this  behavior  at  A).  At  sinkhole  B  only  one  electrode  (No.  21) 
exhibited  a  similar  SP  Increase  over  this  time  period.  These  phenomena  have 
not  been  observed  at  any  of  the  other  sites  during  this  research.  It  is  pos¬ 
sible  that  these  dramatic  but  short-term  events  were  related  to  the  condi¬ 
tions  of  extreme  drought  that  prevailed  at  the  time.  The  Palmer  Index  for  the 
end  of  September  showed  all  regions  of  the  state  in  severe  drought,  according 
to  the  Virginia  Water  News  (Vol  17,  No.  10,  p  6).  Precipitation  at  the 
McHenry  site  totalled  only  13.7  cm  between  21  July  and  13  October.  Figure  17 


depicts  SP  changes  at  sinkhole  A.  The  short-term  anomalies  at  electrodes 
Nos.  26  and  29  can  be  seen  in  the  west  and  southwest  quadrant  from  17  and 
23  September.  On  14  October  there  was  1.3  cm  of  rain,  which  eliminated  these 
two  anomalies.  On  19  October,  after  drying  out,  these  two  electrodes  again 
became  negative,  but  not  to  the  extremes  reached  previously.  These  short-term 
anomalies  occur  in  outer  ring  electrodes  and  apparently  are  not  reflected  by 
corresponding  changes  in  the  inner  ring  of  electrodes.  They  probably  repre¬ 
sent  vertical  flow  or  seepage  of  soil  moisture  into  fractures.  To  sustain 
base  flow  of  ground  water  at  the  spring  (its  output  was  relatively  unaffected 
by  the  drought) ,  water  that  is  tightly  held  in  fractures  and  soil  pores  is 
pulled  downward  into  conduits  and  channels  as  the  more  accessible  storage 
water  is  exhausted  from  the  ground.  This  supply  of  ground  water,  that  is  oth¬ 
erwise  tightly  held  in  storage  by  capillary  forces,  would  likely  generate  high 
streaming  potentials  when  flowing  through  very  small  fractures  (or  sediment- 
filled  fissures).  Perhaps  the  reason  this  effect  was  observed  at  the  McHenry 
sinkholes  is  because  they  are  close  to  base  level  discharge  (the  spring); 
whereas,  all  the  other  sites  are  areas  of  ground-water  recharge.  During 
extreme  drought,  areas  near  a  spring  may  be  subjected  to  conditions  that  make 
extraordinary  demands  on  ground  water  held  in  storage,  and  this  could  create 
ephemeral  anomalies  that  reflect  temporary  areas  of  preferential  flow. 

Whose  Hole  Cave  results 

49.  At  Whose  Hole  Cave  three  different  arrays  of  electrodes  were 
emplaced.  A  ring  of  12  electrodes  was  placed  in  the  sinkhole  containing  the 
cave  entrance  (Table  7) ,  and  three  rows  of  ten  electrodes  were  placed  directly 
above  that  portion  of  the  cave  containing  the  stream  (Table  8).  An  array  of 
32  electrodes  (two  concentric  rings  of  16  each)  was  emplaced  in  December  1986 
(Table  9).  Figure  18  is  a  map  showing  the  360  m  of  passage  that  have  been 
surveyed  and  the  relationship  of  the  electrode  arrays  to  the  underlying  cave 
passages.  Passages  are  strike  controlled  and  trend  to  the  northeast  toward 
the  McHenry  sinkholes  A  and  B.  The  stream  inside  the  cave  has  been  dyed  to 
establish  its  connection  to  the  spring  at  sinkhole  B.  A  weir  and  water  level 
measuring  device  (float  attached  to  a  potentiometer)  provided  streamflow  data 
for  the  cave  stream.  Streamflow  could  not  be  related  to  SP  changes  at  the 
array  of  electrodes  over  the  stream  passage.  This  may  be  due  to  the  depth  of 
25  m  or  to  the  relatively  small  stream  (flows  ranging  from  0.002  m/sec  to 
0.0!  m/sec)  or  to  the  lack  of  sand-sized  sediment  in  the  fill  within  the  cave 
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(mostly  clay  and  gravels).  A  profile  of  average  SP  for  the  array  over  the 
stream  passage  is  shown  in  Figure  19.  One  might  expect  an  anomalously  nega¬ 
tive  area  directly  above  the  cave.  In  this  case  the  reverse  Is  observed — the 
positive  SP  values  being  over  the  cave.  The  two  negative  anomalies  seen  on 
either  side  of  the  array  are  unexpected.  The  anomaly  at  electrodes  Nos.  2  and 
3  may  be  the  result  of  the  intersection  of  a  joint  or  fault  plane  with  the 
surface  here.  Such  a  plane  would  allow  preferential  movement  of  water  that 
would  generate  high  SP  values.  As  can  be  seen  in  Figure  19,  this  plane  is 
exposed  as  a  wall  in  the  cave  passage.  The  cave  contains  numerous  tube-like 
features  that  incline  at  similar  angles  and  which  serve  as  conduits  for 
descending  trickles  of  water.  The  anomaly  at  electrodes  Nos.  7  and  8  cannot 
be  explained.  Perhaps  another  joint  or  fault  that  is  unknown  intersects  the 
surface  here.  The  implications  of  this  data  are  that  SP  results  in  karst 
areas  of  structural  complexity  must  be  interpreted  very  critically  with  regard 
to  underlying  structure.  The  effect  of  inclined  planes  on  preferential  water 
movement  may  mean  that  the  location  of  SP  anomalies  is  offset  from  underlying 
cavities  or  features  of  interest.  This  research  suggests  that  SP  results  over 
Whose  Hole  Cave  are  related  to  soil  water  percolating  downward  into  planes  of 
greater  permeability  and  preferential  flow  rather  than  to  the  flow  of  the 
vadose  cave  stream.  This  result  is  consistent  with  SP  theory  as  expressed  in 
the  literature. 

50.  SP  results  from  the  sinkhole  entrance  to  the  cave  can  be  related  to 
underlying  cave  features.  Plan  views  of  this  portion  of  the  cave  with  the 
associated  electrode  arrays  are  given  in  Figures  20  and  21.  For  Figure  20, 
with  an  array  of  12  electrodes,  the  anomalous  electrodes  are  Nos.  1,  2,  3,  9, 
and  12,  while  all  other  electrodes  have  positive  SP  values.  The  anomalous 
electrodes  are  located  above  cave  passages  where  ground  water  seeps  through 
sediment  fill  (silty  clay  and  gravel  deposits)  or  above  steeply  inclined  nar¬ 
row  fissures  which  channel  downward  percolating  soil  water  into  the  cave. 

These  vertical  fissures  and  tubes  are  relatively  high  in  the  cave  and  extend 
upward  in  proximity  to  the  surface  (especially  electrode  Nos.  1  and  3).  The 
values  for  electrodes  Nos.  9  and  12  are  probably  explained  by  the  fact  that 
the  passage  between  them  contains  a  pool  of  water  during  wet  weather  that  is 
created  by  a  natural  "dam"  of  sediments  and  breakdown.  Water  in  the  pool 
percolates  slowly  through  the  permeable  material  in  the  dam,  moving  from  the 
vicinity  of  electrode  No.  9  toward  electrode  No.  12.  This  direction  of 
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flow  can  be  observed  in  the  cave  but  Is  also  confirmed  by  the  consistently 
greater  SP  values  at  electrode  No.  12  than  at  electrode  No.  9.  During  dry 
weather  there  is  little  or  no  water  feeding  into  this  area  and  the  pool  disap¬ 
pears.  Water  percolating  through  this  dam  of  sediment  would  probably  generate 
large  streaming  potentials.  The  series  of  SP  contour  maps  shown  in  Fig¬ 
ures  22,  23,  and  2A  all  indicate  major  anomalous  zones  within  the  northeast 
quadrant  of  the  array.  Also,  the  contour  map  of  EM  conductivity  (Figure  25) 
shows  the  northeast  quadrant  as  the  principal  EM  anomaly.  The  cave  passages 
under  this  area  are  high,  narrow  fissures  that  are  close  to  the  ground  surface 
serve  as  zones  of  preferential  ground-water  movement.  When  the  EM  conductiv¬ 
ity  survey  was  conducted  in  December  1986,  soil  moisture  content  was  fairly 
high.  It  is  probable  that  the  soil  near  the  tops  of  these  fissure  passages 
(and  the  soil  and  sediment  filling  the  fracture  systems  feeding  these  fis¬ 
sures)  was  saturated  to  a  greater  extent  than  the  surrounding  soil.  There¬ 
fore,  this  zone  of  greater  saturation  would  be  more  conductive  and  would 
produce  the  conductivity  anomaly  seen  in  Figure  25.  Similarly,  the  percola¬ 
tion  of  water  in  this  area  would  generate  anomalous  streaming  potentials  also. 

51.  A  temporal  effect  on  SP  can  be  seen  in  Figure  26.  SP  values  for 
the  12  electrodes  around  the  entrance  sinkhole  are  profiled  for  each  reading 
date,  from  September  through  November.  The  desiccating  effects  of  drought  are 
reflected  as  a  gradual  shift  of  SP  values  from  negative  to  positive.  The  mid¬ 
dle  of  October  showed  all  positive  values  (the  anomalous  spikes  have  even  been 
lost)  which  correspond  to  maximum  drought  conditions.  In  late  October,  as 
some  rain  began  to  return  and  vegetation  became  dormant,  moisture  was 
recharged  to  the  ground,  and  the  SP  profiles  began  to  shift  back  toward  nega¬ 
tive  (with  the  anomalous  peaks  returning).  In  early  November,  conditions  of 
Increasing  seepage  were  observed  in  the  cave  around  electrodes  Nos.  3  and  9. 
The  data  presented  in  Figure  26  appear  to  portray  a  direct  relationship 
between  drought-induced  desiccation  effects,  the  return  of  more  moisture,  and 
SP  values.  Although  electrodes  Nos.  7  and  8  were  located  over  cave  passage, 
they  had  a  consistent  positive  value.  This  was  probably  because  their  pas¬ 
sages  were  relatively  deeper  than  those  areas  beneath  the  anomalous  electrodes 
(about  20  m  as  opposed  to  5  m) .  Also,  even  though  there  is  water  flowing  near 
electrodes  Nos.  7  and  8,  it  is  a  relatively  small  quantity  and  the  flow  is 
over  bare  rock  and  travertine  deposits  (a  condition  that  would  not  be  likely 
to  generate  appreciable  streaming  potentials).  At  this  site  then,  the  SP 
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anomalies  are  In  areas  where  ground  water  percolates  through  a  sediment  dam  or 
plug  and  around  high  vertical  fissures  and  tubes  that  approach  the  surface  and 
serve  as  channels  to  concentrate  the  entry  and  downward  seepage  of  moisture 
from  the  soil  and  eplkarst  Into  the  cave  passages  below.  Passages  about 
15  to  20  m  deep  with  water  flowing  over  rock  or  travertine  do  not  appear  as  SP 
anomalies . 


Time  Varying  Effects  on  SP  Results 

52.  The  time  varying  effects  of  soil  temperature  and  precipitation  on 
SP  values  have  been  observed  at  all  the  study  sites.  It  Is  difficult  to  char¬ 
acterize  the  exact  response  of  a  given  site  or  electrode  to  these  changing 
stimuli  as  the  reactions  sometimes  appear  contradictory  or  erratic  because  of 
an  interactive  effect  of  these  variables.  There  may  be  other  factors  Influ¬ 
encing  SP  values  as  well,  further  complicating  the  detection  of  cause  and 
effect.  Such  factors  would  Include  physical  and  chemical  properties  of  soil, 
depth  to  bedrock,  relief,  vegetation,  aspect  of  the  site,  and  seasonal  cli¬ 
matic  changes.  The  fluctuations  imposed  by  these  factors  constitute  signal 
noise  that  can  mask  or  exaggerate  the  streaming  potential  effects  that  are 
being  sought.  Not  only  do  the  various  sites  respond  differently  to  changing 
environmental  factors,  but  individual  electrodes  within  a  given  array  can 
respond  differently,  as  seen  in  Figure  27.  Here  the  mixed  reactions  of 

four  electrodes  to  3.4  cm  of  rain  are  shown.  It  may  be  that  such  behavior  is 
attributable  to  the  extreme  drought  conditions  that  existed  during  the  study 
period,  and  that  during  periods  of  more  normal  soil  moisture  the  SP  responses 
are  more  predictable. 

53.  Soil  moisture  is  a  critical  parameter  to  consider,  along  with  pre¬ 
cipitation,  as  it  can  be  related  to  soil  saturation  conditions  and  the  lag 
times  between  precipitation  events  and  SP  values.  Soil  moisture  data  also 
compensate  for  the  effects  of  vegetation  types  and  activity  on  infiltration  of 
precipitation.  In  general,  increased  soil  moisture  causes  an  SP  increase 
(values  become  more  negative) ,  presumably  resulting  from  a  greater  streaming 
potential  as  the  water  percolates  through  the  soil  and  regolith.  Figure  28 
illustrates  this  reaction  for  two  test  electrodes.  A  sharp  increase  in  soil 
moisture  produced  an  SP  change  of  grater  than  100  mv  in  the  negative  direc¬ 
tion,  but  this  occurred  about  18  hr  after  the  rain  began.  It  can  also  be  seen 
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that  the  SP  returns  to  its  former  value  very  gradually  over  a  period  of  days. 
However,  depending  on  the  degree  of  saturation,  some  precipitation  events  may 
produce  relatively  minor  SP  changes  or  may  even  result  In  a  positive  change  of 
value . 

54.  The  time  varying  influence  of  temperature  also  changes  SP  values. 
Thermocouple  probes  that  produce  a  digital  reading  on  a  voltmeter  have  been 
used  to  measure  soil  temperature  at  a  depth  of  about  15  cm.  The  tip  of  the 
electrodes  are  deep  enough  to  be  at  a  constant  or  slowly  changing  temperature; 
whereas,  the  exposed  and  near-surface  area  is  subject  to  fluctuations  of  up  to 
10°  C.  A  diurnal  cycle  in  the  SP  values  has  been  observed  that  corresponds  to 
the  warming  and  cooling  effect  of  the  surface  soil.  This  effect  is  more  pro¬ 
nounced  when  soil  moisture  is  at  low  levels,  as  can  be  seen  in  Figure  28. 
Similarly,  greater  daily  temperature  extremes  produce  a  larger  SP  cyclic 
change.  If  the  exposed  portion  of  the  electrode  is  warmed  by  the  sun,  partic¬ 
ularly  on  cold  mornings,  a  large  SP  change  can  result  from  the  thermocouple 
effect.  Such  changes  are  systematic  and  affect  the  SP  values  relatively  uni¬ 
formly  for  a  site.  An  SP  change  of  50  mv  has  been  observed  for  frosty  morn¬ 
ings  followed  by  sunny  days  that  produced  a  surface  soil  temperature  range  of 
4°  C  and  an  air  temperature  range  of  16°  C.  Temperature  effects  on  SP  are  not 
yet  fully  understood,  as  evidenced  by  the  observation  that  at  some  sites  SP 
increases  as  temperature  increases,  and  at  other  sites  SP  decreases  as  temper¬ 
ature  Increases.  More  work  needs  to  be  done  to  understand  the  complex  inter¬ 
action  between  temperature,  moisture,  SP,  and  other  factors  which  may  be  site 
specific. 

55.  Until  such  effects  can  be  isolated  and  predictably  quantified, 
using  SP  results  based  on  averages  of  many  readings  may  be  the  only  practical 
approach  for  describing  the  general  response  of  a  site  to  variables.  To 
understand  the  nature  of  individual  electrode  responses  to  processes  affecting 
it,  long-term  accurate  monitoring  is  required.  At  the  long-term  monitoring 
site,  data  were  collected  from  1  August  1986  to  8  October  1986  using  volt¬ 
meters  and  thermometers  (Tables  10,  11,  12,  and  13).  Data  at  this  site  were 
collected  automatically  by  the  ADC  system  from  8  October  1986  to  11  May  1987 
(Tables  14-19).  The  automatic  data  collection  system  described  previously  has 
great  potential  for  such  application.  Temperature  and  moisture  sensors  can  be 
placed  in  the  ground  and  data  can  be  sampled  at  rates  adequate  for  detecting 
subtle  changes.  The  data  presented  in  Figures  28  and  29  were  recorded  with 
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the  ADC  at  the  electrode  monitoring  site  (dry  stream  bed  below  Whose  Hole 
Cave).  From  Figure  28,  the  data  show  soil  moisture  responding  to  diurnal 
temperature  cycles,  but  the  SP  record  does  not  show  a  similar  response  when 
moisture  levels  are  high.  Also,  when  soil  moisture  Is  high,  the  SP  changes, 
as  a  result  of  additional  precipitation,  are  lessened  by  a  factor  of  about  two 
(60-mv  change  on  5  November  as  opposed  to  140-mv  change  on  25  October)  com¬ 
pared  with  dry  soil  conditions. 

56.  During  the  eight  month  period  that  the  ADC  system  was  used  to  mea¬ 
sure  SP  at  the  five  electrodes  of  the  long-term  monitoring  site,  electrode 
No.  5  was  the  only  electrode  that  was  always  negative.  The  location  of  elec¬ 
trode  No.  5  did  not  seem  to  be  any  different  from  the  other  electrodes  except 
that  the  soil  there  appeared  to  have  more  clay.  Consequently,  soil  samples 
were  taken  by  coring  beside  each  of  the  five  electrodes,  and  a  partlcle-slzed 
distribution  was  determined  by  sieve  and  hydrometer  analysis.  The  resulting 
gradation  curves  are  shown  in  Figure  30  and  It  can  be  seen  that  the  soil  at 
electrode  No.  5  does  Indeed  contain  appreciably  more  clay  than  the  other  elec¬ 
trodes.  The  coarsest  material  is  the  gravelly  soil  from  the  streambed  at 
electrode  No.  3  (which  was  positive  when  dry  and  negative  when  the  Intermit¬ 
tent  stream  flowed) .  The  relationship  between  high  clay  content  and  negative 
SP,  as  shown  by  electrode  No.  5,  is  of  Interest  because  this  seems  to  contra¬ 
dict  the  SP  literature,  which  reports  positive  charges  associated  with  clay 
and  negative  streaming  potentials  associated  with  sand-sized  particles. 
Although  this  is  only  one  data  point  represented  here.  It  seems  to  be  sup¬ 
ported  by  subjective  observations  In  the  field.  While  pulling  electrodes  from 
the  ground  at  various  sites,  many  electrodes  that  had  been  consistently  anoma¬ 
lous  (negative)  were  very  difficult  to  extract  because  of  their  being  in  a 
sticky,  plastic  clay.  It  may  be  that  the  greatest  negative  anomalies  in  karst 
terrain  are  not  a  result  of  streaming  potentials  in  sand-sized  material  but. 
Instead,  are  from  Ionic  charges  or  membrane  activity  in  clay  minerals.  The 
mineralogy  of  carbonate  residual  soils  could  have  a  drastically  different  SP 
effect  than  that  of  soils  derived  from  siliceous  lithologies,  and  would  there¬ 
fore  impose  different  criteria  for  interpreting  ground-water  flow  regimes  from 
SP  data.  Also,  the  surflcial  aspect  of  a  short  electrode  stuck  into  a  high 
clay  content  soil  ought  to  be  considered.  If  the  clay  overlaid  a  sandy  sub¬ 
soil  or  permeable  regolith  through  which  ground  water  was  flowing,  it  could  be 
possible  that  the  ionic  charges  associated  with  the  clay  could  mask  or  amplify 
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any  streaming  potential  charges  produced  at  a  greater  depth.  The  vertical 
aspect  of  SP  phenomena  is  an  area  that  needs  to  be  explored,  along  with  the 
other  factors  previously  mentioned. 


Sinkhole  Artificial  Recharge  Investigations 

57.  Since  the  Moore  sinkhole  A  was  an  unsuccessful  pond  site,  it  pre¬ 
sumably  contains  one  or  more  leakages  and  would  therefore  be  an  ideal  place  to 
conduct  an  experiment  in  detecting  flow  paths  and  rates.  An  array  of  32  elec¬ 
trodes  (two  concentric  rings  with  diameter  of  15  and  45  m)  was  installed.  SP 
data  were  collected  by  voltmeter  (Table  20).  Two  separate  tests  were  then 
conducted.  For  each  test,  water  trucks  discharged  a  total  of  15,900  I  into 
the  sinkhole  over  a  period  of  4  to  5  min.  SP  was  monitored  before,  during, 
and  after  the  discharge,  both  by  the  ADC  system  and  manually  by  voltmeters. 

58.  Both  the  water-dumping  tests  were  successful  in  demonstrating  that 
a  negative  SP  shift  occurred  in  response  to  the  water  draining  through  the 
sinkhole  and  that  the  greatest  SP  anomaly  occurred  around  the  same  electrode 
(No.  3)  on  both  occasions.  Although  moisture  and  temperature  conditions  dif¬ 
fered  for  the  two  tests,  the  No.  3  electrode  produced  a  strong  anomaly  in  SP 
both  times,  as  shown  in  Figure  31.  Since  electrode  No.  3  is  located  on  the 
inner  ring  of  electrodes,  one  could  expect  the  opposite  electrode  in  the  outer 
ring  (No.  19)  to  also  show  a  peak  if  the  flow  is  outward  in  a  lateral  plane. 
Since  this  is  not  the  case,  it  appears  that  the  water  was  descending  verti¬ 
cally  in  the  vicinity  of  No.  3,  and  was  not  reaching  the  outer  ring  of 
electrodes. 

59.  For  the  first  test  on  4  October  1986,  the  ground  was  extremely  dry 

from  drought  conditions  (the  preceding  2  months  had  received  only  8  cm  of 

rain,  with  no  rain  at  all  for  the  month  before  the  test) .  Diurnal  temperature 

swings  were  not  great  at  that  time.  The  discharged  water  ponded  in  the  bottom 

of  the  sinkhole  to  a  depth  of  about  10  cm  and  took  70  min  to  be  absorbed. 

This  represents  a  percolation  rate  of  roughly  227  i/mln  over  an  area  of  about 
2 

177  m  .  The  Induced  SP  changes  are  shown  as  a  polar  plot  in  Figure  32.  Elec¬ 
trode  No.  3  experienced  an  SP  change  of  about  200  mv  and  goes  off  the  scale 
compared  with  the  relatively  smaller  changes  of  the  other  electrodes  (ranging 
from  0  to  32  mv) .  The  time  required  to  reach  the  maximum  change  varied  from 
2  to  24  hr  for  the  various  electrodes. 
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60.  During  the  first  test,  the  water  was  dumped  in  the  northeast  quad¬ 
rant  of  the  sinkhole  and  during  the  discharge  process  water  flowed  over  elec¬ 
trodes  Nos.  1,  2,  and  3.  Their  response  was  immediate,  as  can  be  seen  by  the 
charts  in  Figure  33.  This  illustrates  the  initial  response  to  wetting  which 
occurs  in  a  matter  of  seconds.  The  sudden  sharp  negative  increase  followed  by 
a  sharp  positive  shift  is  probably  the  "moto-electric"  effect  that  is 
described  by  Ogllvy,  Ayed,  and  Bogoslovsky  (1969).  This  effect  arises  when  a 
metallic  electrode  is  washed  by  the  streaming  potential  fluid  (water,  in  this 
case),  which  causes  a  negative  SP  shift,  the  magnitude  of  which  is  a  function 
of  the  metal.  Copper  is  given  a  value  of  -26  mv,  which  is  about  the  same 
value  as  the  recovery  period  (positive  shift)  seen  for  electrode  No.  3.  Like¬ 
wise,  electrode  No.  2  went  from  an  initial  value  of  -8  mv  to  one  of  +25  mv  (a 
total  of  33  mv  change  due  to  moto-electric  effect)  before  settling  into  a 
steady  negative  Increase. 

61.  The  second  test  was  conducted  2  weeks  later  on  18  October  1986. 

The  soil  was  more  moist  than  during  the  first  test,  since  1.8  cm  of  rain  had 
fallen  4  days  before.  Temperature  extremes  were  much  greater  because  of  heavy 
frost  followed  by  clear  sunny  days.  Figure  34  illustrates  the  temperature 
ranges  of  16®  C  for  air  and  4®  C  for  surface  soil.  The  corresponding  plots  of 
SP  values  show  a  pronounced  diurnal  cycle  because  of  the  influence  of  tempera¬ 
ture.  This  effect  produced  a  change  in  SP  of  more  than  50  mv  in  a  systematic 
fashion  that  served  to  mask  the  subtle  SP  changes  that  may  have  been  induced 
by  the  water  discharge,  with  the  exception  of  electrode  No.  3.  Preferential 
water  movement  at  electrode  No.  3  was  again  Implied  by  the  anomalous  SP  change 
of  94  mv,  nearly  twice  the  value  of  all  the  other  electrodes.  A  polar  plot  of 
the  change  in  SP  for  the  second  test  is  nearly  circular  with  the  exception  of 
a  peak  at  No.  3.  For  the  second  test,  the  water  was  discharged  from  the 
northwest  side  of  the  sinkhole  so  that  electrode  No.  3  would  not  be  directly 
wetted  as  occurred  during  the  first  test.  Data  were  collected  by  the  ADC  sys¬ 
tem  for  a  day  following  the  test,  and,  as  seen  in  Figure  34,  it  discloses  the 
presence  of  the  diurnal  temperature  cycle  influencing  the  SP  values  in  a 
nearly  identical  manner  for  both  days.  The  anomalous  behavior  of  electrode 
No.  3  to  the  water  discharge  tests  is  also  evident  from  the  SP  record  shown  in 
Figure  35.  Both  tests  produced  significant  negative  shifts  at  No.  3,  but  only 
produced  a  relatively  minor  negative  shift  for  the  array  average  SP. 
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62.  During  the  second  test  it  took  165  min  for  the  15,900  I  of  water  to 
percolate  down  into  the  soil,  a  rate  of  95  £/min  for  the  same  area  as  in 

test  1.  This  slower  rate  was  presumably  because  of  the  greater  moisture  con¬ 
tent  of  the  soil  during  the  second  test.  Although  flow  rates  for  the  second 
test  are  less  than  half  the  value  during  the  first  test,  it  is  difficult  to 
equate  the  figures  to  SP  changes  because  of  the  large  systematic  Influence  of 
temperature  that  masked  most  of  the  SP  data.  However,  it  Is  interesting  to 
note  the  correspondence  between  SP  change  of  electrode  No.  3  and  the  flow  rate 
estimates — a  ratio  of  about  2:1  in  both  cases: 

Flow  Rates  ASP  (mv) 

(Z/mln)  No.  3 

Test  1  227  19A 

Test  2  95  94 

63.  From  Figure  32  it  can  be  seen  that  the  location  of  electrode  No,  3 
is  in  the  northeastern  quadrant  of  the  sinkhole.  This  is  on  the  down-gradient 
side  of  the  sinkhole,  and  it  is  reasonable  to  assume  that  water  is  probably 
leaving  the  sinkhole  in  this  direction,  which  corresponds  to  the  strike  of  the 
underlying  beds.  Figure  32  is  a  polar  plot  of  azimuthal  electrical  resistiv¬ 
ity  data  on  several  dates,  showing  an  ellipse  oriented  to  the  northeast  also. 
This  data  suggests  that  resistivity  of  the  subsurface  material  is  greater 
along  the  direction  of  strike  than  it  is  at  right  angles  to  the  strike.  A 
possible  explanation  for  this  phenomenon  may  be  that  proportionally  more  fis¬ 
sures  or  conduits  are  developed  in  the  strike  direction  (most  cave  passages  in 
the  area  are  strike  oriented) ,  and  these  are  mainly  empty  so  that  the  collec¬ 
tive  volume  of  air-filled  spaces  increases  the  resistivity.  If  this  model  is 
correct,  the  resistivity  along  the  northeast  axis  should  decrease  as  more 
water  enters  and  fills  these  channels  (water  being  more  conductive  than  air, 
of  course).  This  response  does  take  place  for  the  two  dates  plotted.  The 

18  October  data  were  collected  several  hours  after  the  water  was  discharged 
into  the  sinkhole.  The  data  for  24  October,  nearly  a  week  (with  no  rain) 
later,  show  that  resistivity  has  Increased  along  the  northeast  axis  of  the 
ellipse,  suggesting  that  some  water  has  drained  out  of  the  fissures  and 
fractures. 

64.  Following  the  two  water  discharge  tests  in  the  Moore  sinkhole  A, 
the  investigation  at  this  site  was  expanded  to  include  two  other  sinkholes 
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(B  and  C)  which  are  in  close  proximity  (Table  21).  The  general  vicinity  is  a 
complex  of  sinkholes,  and  it  was  hoped  that  water  draining  out  of  sinkhole  A 
(the  highest  in  elevation  of  the  three)  would  be  detectable  as  a  flow  path 
leading  to  either  or  both  of  the  sinkholes  B  and  C.  In  addition  to  SP  mea¬ 
surements,  an  EM  conductivity  survey  was  made  at  the  three  sinkholes  and  the 
results  are  shown  in  Figures  36-39.  Sinkholes  A  and  C  appear  to  be  more 
closely  related  to  one  another  in  terms  of  conductivity  anomalies.  As  seen  in 
Figure  39,  the  anomalous  zones  of  A  and  C  are  also  related  to  topographic 
expression  as  well,  being  clustered  at  the  point  where  surface  water  flow 
would  enter  sinkhole  A  and  seep  downward  in  the  center  of  the  sinkhole,  then 
again  where  surface  water  would  first  concentrate  in  the  western  limb  of  sink¬ 
hole  C  and  flow  subterraneously  out  of  C  to  the  northeast  toward  another  large 
sinkhole.  This  northeasterly  flow  direction  out  of  C  also  corresponds  to  the 
geologic  strike.  At  sinkhole  B  the  conductivity  anomaly  is  located  in  the 
northwest  quadrant  where  a  small  (2  m)  sinkhole  occurs  on  the  edge  of  B.  This 
anomaly  also  aligns  on  the  northeast  (strike)  direction  from  the  nearby  pond. 
Based  on  the  conductivity  data,  one  could  Interpret  a  ground-water  drainage  or 
flow  path  entering  sinkhole  A  on  the  west,  descending  in  the  center  of  the 
sinkhole  and  moving  eastward  to  sinkhole  C.  In  sinkhole  C  the  flow  then 
appears  to  be  directed  along  strike  toward  the  next  sinkhole  adjacent  to  the 
northeast.  The  locational  relationships  of  the  three  Moore  sinkholes  are  seen 
in  Figure  40. 

65.  The  SP  data  from  the  Moore  sinkhole  appear  to  be  compatible  with 
the  interpretation  based  on  EM  conductivity.  Figures  41,  42,  and  43  show  SP 
contours  at  sinkholes  A,  B,  and  C,  respectively.  These  SP  data  were  collected 
in  early  January  1987,  a  few  days  after  the  EM  conductivity  survey  was  made  in 
later  December  1986.  The  SP  anomalies  largely  agree  with  the  conductivity 
anomalies  so  that  the  same  interpretation  of  a  flow  path  from  A  to  C  seem  rea¬ 
sonable.  In  sinkhole  B  (Figure  42)  a  very  obvious  SP  anomaly  occurs  above  the 
small  2-m  sinkhole  in  the  same  location  as  the  EM  conductivity  anomaly  of  Fig¬ 
ure  37.  Since  Moore  sinkhole  A  has  been  monitored  for  SP  for  several  months 
before  sinkholes  B  and  C  were  studied,  these  data  are  shown  in  Figure  44.  The 
SP  anomalies  are  slightly  different  for  this  time  period  (which  occurred  dur¬ 
ing  the  drought) ,  but  they  still  tend  to  suggest  a  flow  direction  from  west  to 
east  and  toward  sinkhole  C.  All  of  this  data  show  a  relationship  between  SP 
and  EM  conductivity,  and  it  also  suggests  a  ground-water  flow  connection 
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between  sinkholes  A  and  C.  What  It  does  not  help  explain  Is  the  anomalous  SP 
response  of  electrode  No.  3  In  sinkhole  A  during  the  two  water  discharge 
tests.  This  electrode.  In  the  northeast  quadrant  of  sinkhole  A,  does  not 
appear  anomalous  In  any  of  the  other  SP  or  EM  conductivity  data.  Perhaps  a 
combination  of  the  drought  conditions  that  prevailed  In  October  1986  and  the 
fact  that  flow  was  artificially  Induced  offer  some  explanation  for  the 
response  of  electrode  No.  3  at  that  time.  More  water  discharge  tests,  exca¬ 
vating  sinkhole  A  to  bedrock,  or  penetrometer  tests  around  electrode  No.  3  may 
be  the  only  way  to  ascertain  what  Is  really  happening.  Despite  this  contra¬ 
diction,  the  bulk  of  the  data  from  the  Moore  site  suggests  that  SP  anomalies 
around  sinkholes  can  be  used  to  Infer  flow  paths,  and  that  EM  conductivity 
seems  to  confirm  this.  These  Inferred  flow  paths  are  probably  shallow — ^within 
a  20-m  depth,  as  the  EM  data  at  10  and  20  m  agree  with  the  Interpretation,  but 
the  40-m  survey  appears  rather  amorphous  In  this  regard. 

66.  Other  results  from  the  long-term  monitoring  of  SP  electrodes  and 
environmental  variables  are  shown  In  Figures  45-56  and  In  Tables  13-22.  Ini¬ 
tially  these  data  were  collected  manually  (Tables  13,  14,  15)  from  1  August 
1986  through  8  October  1986,  after  which  the  automated  data  collection  (ADC) 
system  was  employed  until  termination  on  11  May  1987.  The  graphs  shown  In 
Figures  45-56  were  generated  by  computer  (Lotus  program)  using  the  digital 
data  collected  automatically  and  do  not  reflect  the  time  period  1  August  - 

8  October.  Of  the  five  SP  electrodes  that  were  monitored.  No.  5  was  always 
negative.  No.  4  was  always  positive.  No.  2  began  positive  but  after  November 
It  changed  to  negative  for  the  remaining  monitored  time  period,  and  No.  1 
varied  between  negative  and  positive  values  over  the  entire  time  period. 
Electrode  No.  3  (Figure  47)  was  positive  as  long  as  the  Intermittent  streambed 
In  which  It  was  located  remained  dry,  but  when  the  stream  flowed  temporarily 
after  heavy  rainfall  the  SP  values  became  negative,  as  seen  in  the  sharp  peaks 
in  Figure  52.  Presumably,  these  anomalous  SP  values  are  a  result  of  streaming 
potentials  generated  by  the  flow  of  water  through  the  sediments  of  the  stream 
bed.  Figure  55  shows  standard  deviation  of  electrode  No.  3  and  the  greatest 
values  here  are  because  of  the  extreme  changes  in  SP  (within  a  24-hr  period) 
as  a  result  of  this  streaming  potential. 

67.  The  five  electrodes  at  the  long-term  monitoring  site  were  monitored 
at  sampling  rates  of  one-half  hour  to  4  hr,  and  these  data  have  been  averaged 
over  24  hr  to  produce  the  data  In  Tables  13-22  and  Figures  45-56.  Changes  In 
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standard  deviation  shown  are  thus  relative  to  SP  variations  within  a  24-hr 
period.  Standard  deviation  and  SP  change  are  proportional  as  seen  In  Fig¬ 
ure  54.  All  of  the  five  electrodes  show  similar  SP  responses  initially 
(October),  but  the  first  significant  rainfall  to  occur  after  implanting  the 
electrodes  (1.35  In.  on  25  October)  caused  a  major  negative  shift  in  all  of 
them,  after  which  their  SP  variations  were  dissimilar.  Throughout  the  8-month 
period,  episodes  of  major  rainfall  produced  marked  negative  SP  responses  In 
all  five  electrodes. 

68.  The  effect  of  soil  moisture  (Figure  51)  on  SP  response  is  less 
obvious  than  that  of  rainfall.  Soil  moisture  was  monitored  electronically 
(the  scale  in  Figure  51  Is  in  relative  "digital"  units  where  high  values 
represent  higher  moisture  content)  from  a  gypsum  block  sensor  buried  at  20-cm 
depth  in  the  soil.  The  pattern  of  soil  moisture  in  Figure  51  is  similar  to 
that  of  soil  temperature  (Figure  50),  but  this  may  be  from  a  thermal  effect 
transmitted  from  wires  exposed  to  atmospheric  temperature  changes.  There  are 
undoutably  lag  times  between  rainfall  and  soil  molsture/SP  effects,  but  these 
are  not  immediately  obvious  in  the  data.  The  extreme  fluctuation  of  soil 
moisture  data  shown  In  the  beginning  of  Figure  51  Is  because  of  artificial 
effects  produced  by  dumping  water  on  the  sensor  to  test  Its  performance. 

69.  Soil  temperature  seems  to  play  a  more  direct  role  In  affecting  SP 
than  does  soil  moisture.  Taken  alone,  soil  temperature  seems  to  be  Inversely 
proportional  to  SP  (at  this  site),  as  seen  in  Figure  53.  In  this  figure  the 
SP  record  for  electrode  No.  2  has  been  inverted  and  compared  with  the  soil 
temperature  record.  They  agree  favorably  over  most  of  the  time  period,  but 
the  pattern  Is  reversed  under  extremely  wet  conditions  as  seen  in  late  March 
through  April.  This  Implies  that  soil  temperature  effects  are  masked  or 
reversed  when  soil  moisture  reaches  a  threshold  level.  This  emphasizes  the 
complex  Interactive  nature  of  environmental  variables  on  SP.  Extreme  condi¬ 
tions  In  the  environmental  variables  will  override  the  other  factors  and  domi¬ 
nate  or  reverse  SP  responses. 
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PART  IV:  SUMMARY  AND  CONCLUSIONS 


70.  Spontaneous  potential  has  been  monitored  at  a  variety  of  karst  fea¬ 
tures  for  a  period  of  5-9  months.  Significant  SP  anomalies  were  measured  at 
all  sites,  and  are  assumed  to  be  due  primarily  to  streaming  potentials  created 
by  ground-water  flow.  The  magnitude  of  SP  anomalies  and  SP  changes  is  related 
to  ground-water  flow  rates  along  with  a  host  of  other  factors  which  make 
interpretation  difficult.  Direction  of  flow  is  Indicated  by  concentration  of 
SP  gradients.  Flow  paths  are  inferred  by  SP  data  at  all  of  the  sites  and  seem 
especially  well  defined  at  the  Moore,  Spring  Branch,  McHenry,  and  Harris 
sites.  Geologic  studies  along  with  electrical  resistivity  and  EM  conductivity 
surveys  support  the  interpretation  of  flow  paths  at  these  sites. 

71.  The  detected  SP  anomalies  at  these  sites  appear  to  be  due  primarily 
to  water  seepage  and  percolation  through  porous  material  rather  than  conduit 
or  fracture  flow.  There  is  some  evidence  from  Spring  Branch  for  a  general 
systematic  registering  of  conduit  flow  from  the  electrode  array,  but  individ¬ 
ual  electrodes  could  not  be  related  to  changing  flow  regimes.  The  stream  in 
Whose  Hole  Cave  could  not  be  detected  by  the  electrodes  overhead.  SP  probably 
provides  indirect  evidence  for  conduit  flow  by  indicating  anomalies  produced 
by  the  seepage  of  water  from  soil  and  regolith  into  small  fractures  that  sup¬ 
ply  water  to  nearby  conduits.  It  seems  possible  that  SP  can  distinguish 
between  lateral  and  vertical  flow,  as  seen  at  the  Ford  and  Moore  sites.  SP 
was  able  to  detect  underlying  cave  passages  at  Whose  Hole  Cave  if  they  were 
relatively  shallow  (5-10  m)  and  if  they  contained  water  seeping  through  sedi¬ 
ment  fill  or  else  intercepted  water  seepage  from  the  soil  and  regolith 
storage . 

72.  On  the  basis  of  water  discharge  tests  at  Moore  sinkhole  there  seems 
to  be  promise  of  relating  flow  rates  to  flow- induced  SP  changes.  At  the  sites 
where  conduit  flow  is  known  to  exist  (Spring  Branch,  Fords,  Whose  Hole  Cave), 
it  has  not  been  possible  to  correlate  stream  flow  with  SP. 

73.  The  relationship  between  SP  and  environmental  variables  is  complex 
and  not  well  understood.  Soil  temperature,  soil  moisture,  and  precipitation 
are  the  primary  factors  that  have  been  examined,  although  this  research 
Implies  that  the  Influence  of  soil  properties  such  as  clay  content  or  mineral¬ 
ogy  may  also  be  a  dominant  factor  in  SP,  at  least  in  karst  areas.  The  inter¬ 
active  nature  of  these  Influences,  and  perhaps  others,  requires  accurate 
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long-term  monitoring  to  isolate  and  characterize  them.  If  causal  relation¬ 
ships  and  systematic  variations  can  be  identified  and  quantified,  it  would 
allow  more  precision  in  interpreting  and  applying  SP  results.  ADC  systems 
employing  computers  and  remote  sensors  are  a  means  of  gathering  sufficient 
data  to  attack  the  problems  presented  by  environmental  variables.  Such  a  sys¬ 
tem  was  Implemented  in  the  course  of  this  research  and  has  functioned  success¬ 
fully.  Refinements  of  this  system  will  allow  useful  applications  to  long-term 
seepage/flow  monitoring.  Of  theoretical  benefit,  the  nature  of  recharge  and 
storage  processes  in  karst  aquifers  could  be  better  explained  with  the  moni¬ 
toring  capabilities  inherent  in  an  ADC  system.  Some  practical  applications 
are  geotechnical  in  nature;  for  example,  foundation  studies  and  stabilization 
techniques  can  be  monitored  and  supplemented,  well  pump  tests  can  indicate 
aquifer  characteristics  from  SP  effects,  and  long-term  SP  monitoring  can 
indicate  development  or  Increase  of  anomalous  seepage  from  water  retention 
structures.  Possible  environmental  applications  may  exist  in  areas  of  salt¬ 
water  intrusion  and  contamination  plumes  of  high  conductivity  emanating  from 
waste  sites. 

74.  In  summary,  the  SP  technique  is  an  effective  method  for  detecting 
and  mapping  shallow  flow  paths  in  soil,  regollth,  and  rock  around  most  karst 
features.  Its  utility  in  detecting  deep  flow  through  conduits  is  not  resolved 
in  the  current  study.  SP  can  possibly  be  related  to  flow  rates  if  the  flow  is 
primarily  through  porous  material.  Computer  technology  promises  to  further 
develop  the  SP  method  for  successful  applications  to  hydrogeologic  problems  in 
karst  areas  and  is  addressed  in  Reports  3  and  4  in  this  series. 
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Table  1 


Data  from  ADC  System  for  SP  Electrode  Monitoring  Site 


November  4-6,  1986  (Data  used  to  construct  Figure  28 
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ADC  Data  of  SP  and  Soil  Temperature  at  Moore  Sinkhole  During  Water  Dumi 
Test  on  October  4,  1986  (Data  used  to  develop  Figures  31  and  32) 
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Table  3 

SP  Values  (millivolts)  at  the  Spring  Branch  Sinking  Stream  Site 
(Data  collected  by  voltmeter  between  July  1986  and  April  1987) 
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SP  Values  (millivolts)  at  McHenry  Sinkhole  A  from  July  1986  to  April  1987 
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Table  6 

SP  Values  (millivolts)  at  the  McHenry  Sinkhole  B  for  July  1986  to  April  1987 
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Table  7 

SP  Values  (In  millivolts)  for  12  Electrode  Array  at  Cave  Sink 
for  September  to  December  1986 
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Table  9 

'  SP  Values  (millivolts)  for  the  32  Electrode  Array  at  Cave  Sink 

>  (Data  collected  by  voltmeter  between  December  1986 

>  and  April  1987) 


DATE 

12/21 

12/25 

1/13 

2/5 

3/3 

4/28 

SOIL  TEMP  (C) 

14 

4 

2 

6 

9 

FLOW  (cfs) 

.  17 

.  86 

• 

.  5 

.  9 

ELECTRODE  # 

1 

-15 

66 

“*5 

30 

46 

n 

-200 

-168 

-289 

—292 

-294 

-277 

3 

—  12'^ 

-155 

-123 

-152 

-174 

-182 

4 

-175 

-139 

-235 

-265 

-247 

-318 

5 

-156 

-128 

-233 

-230 

-240 

-266 

6 

-302 

oo8 

-313 

-268 

-297 

7 

55 

100 

50 

57 

50 

219 

8 

-265 

-181 

-250 

-266 

-252 

-298 

9 

118 

201 

127 

90 

63 

99 

10 

197 

255 

290 

321 

336 

341 

1 1 

-74 

-14 

-85 

-82 

—  1  ‘J:*3 

-177 

12 

170 

177 

103 

114 

112 

49 

13 

217 

234 

177 

160 

185 

169 

14 

204 

229 

157 

161 

189 

204 

15 

121 

121 

27 

.ji8 

41 

36 

16 

-70 

-60 

-99 

-98 

-113 

-97 

17 

1 10 

137 

100 

130 

1 2'^ 

100 

19 

-10 

-17 

-86 

-44 

-45 

- 1 20 

19 

-300 

-277 

-303 

-294 

—29'^ 

-274 

20 

-240 

-166 

-183 

-194 

-167 

-188 

21 

-280 

-247 

-321 

-304 

-259 

-264 

22 

63 

118 

15 

44 

94 

328 

23 

81 

134 

171 

256 

245 

297 

24 

-32 

-11 

_ 'TT 

-42 

-10 

-46 

25 

173 

203 

199 

217 

233 

”729 

26 

90 

149 

175 

200 

203 

214 

27 

35 

84 

50 

15 

35 

21 

28 

150 

191 

216 

230 

266 

298 

29 

275 

320 

311 

324 

337 

337 

30 

226 

254 

258 

276 

284 

276 

31 

194 

149 

7 

62 

68 

62 

-TO 

198 

198 

163 

91 

85 

68 

I 

I 


Table  10 

SP  Values  in  Millivolts  for  Long-Term  Monltorln; 
Array  of  Electrodes 


ELECTRODES 

SOIL  TEMP 

DATE 

TIME 

#1 

#2 

#3 

#4 

#5 

6" 

1 

AUG 

1100 

251 

121 

233 

187 

214 

23 

1200 

269 

120 

269 

210 

192 

24 

1300 

290 

129 

302 

235 

172 

24 

1 400 

304 

140 

323 

252 

171 

24 

1500 

313 

150 

338 

264 

167 

24 

1700 

327 

167 

358 

284 

169 

24 

1800 

33 1 

171 

368 

295 

164 

25 

1900 

334 

176 

376 

300 

162 

24 

2000 

334 

184 

378 

302 

159 

24 

2300 

314 

182 

375 

292 

127 

23 

2 

AUG 

0700 

316 

198 

399 

314 

101 

21 

0800 

328 

209 

408 

328 

100 

21 

1 000 

343 

224 

433 

349 

115 

23 

1200 

349 

442 

358 

117 

24 

1400 

355 

243 

444 

368 

120 

24 

1 600 

364 

260 

454 

•j83 

124 

24 

1800 

367 

267 

456 

389 

1 22 

24 

2000 

369 

273 

460 

393 

121 

24 

3 

AUG 

0700 

341 

263 

446 

386 

85 

21 

0900 

340 

262 

448 

307 

05 

22 

1 100 

353 

275 

464 

404 

92 

23 

1 300 

361 

285 

471 

416 

102 

24 

1500 

365 

287 

472 

419 

101 

25 

1700 

366 

292 

478 

422 

102 

25 

1900 

367 

293 

457 

424 

100 

25 

2100 

343 

280 

446 

402 

82 

23 

4 

AUG 

0600 

320 

259 

436 

391 

49 

21 

0800 

337 

272 

457 

41 1 

62 

21 

1 000 

348 

280 

470 

425 

74 

22 

1500 

358 

289 

473 

435 

81 

23 

1700 

359 

289 

473 

437 

80 

23 

1900 

358 

289 

474 

436 

78 

23 

2100 

352 

288 

465 

433 

73 

23 

5 

AUG 

0600 

324 

264 

447 

415 

43 

22 

0800 

333 

271 

457 

426 

52 

22 

1700 

356 

290 

400 

450 

72 

23 

2200 

345 

287 

466 

442 

66 

23 

6 

AUG 

0600 

326 

269 

450 

432 

41 

22 

1 700 

60 

55 

1 19 

127 

4 

23 

1 900 

65 

61 

115 

1  oo 

-6 

23 

2 1 00 

66 

50 

1 12 

140 

-20 

22 

2300 

67 

69 

112 

145 

22 

7 

AUG 

0600 

66 

72 

119 

156 

-26 

21 

0800 

65 

71 

120 

155 

-27 

21 

2000 

74 

80 

147 

164 

-14 

23 

8 

AUG 

<;)700 

66 

76 

163 

171 

-31 

23 

1000 

67 

77 

166 

173 

-29 

'7'? 

1300 

69 

78 

168 

174 

-24 

23 

1 800 

73 

92 

173 

179 

-16 

RAINFALL 


.B"  @1500 


Table  11 

SP  Values  in  Millivolts  for  Long-Term  Monltorln; 
Array  of  Electrodes 


ELECTRODES  SOIL  TEMP  RAINFALL 


DATE 

TIME 

#1 

#2 

#3 

#4 

#5 

6" 

9  AUS 

0800 

67 

75 

177 

182 

-40 

21 

1000 

69 

75 

179 

184 

—35 

21 

1200 

77 

83 

188 

192 

-20 

23 

2000 

82 

88 

194 

198 

-19 

24 

10  AUG 

0800 

81 

88 

202 

205 

—33 

or* 

1300 

91 

96 

211 

214 

-12 

2S 

2000 

93 

98 

21 1 

215 

-17 

24 

0.  1" 

1900 

102 

109 

230 

233 

-14 

24 

12  AUG 

0800 

104 

112 

239 

241 

-23 

22 

1900 

128 

1 36 

265 

268 

1 

23 

13  AUG 

0700 

122 

132 

263 

266 

-16 

21 

1900 

147 

157 

291 

293 

7 

22 

14  AUG 

0800 

143 

154 

291 

292 

21 

18  AUG 

1600 

141 

166 

276 

289 

4 

23 

.5" (18  AUG) 

20  AUG 

1500 

-34 

-41 

70 

58 

-263 

21 

.6" ( 19  AUG) 

1800 

-18 

-41 

63 

42 

-260 

21 

1.1" (20  AUG) 

2000 

-12 

-47 

65 

44 

-257 

21 

.  4" 

21  AUG 

0800 

30 

-84 

66 

8 

-262 

21 

.2" 

1800 

51 

-96 

82 

—2 

-254 

22 

22  AUG 

0900 

77 

-132 

87 

-43 

-259 

21 

1800 

85 

-121 

98 

—36 

—242 

24 

23  AUG 

0800 

79 

-132 

114 

-32 

-210 

21 

1800 

81 

-109 

118 

-8 

-193 

25 

24  AUG 

0900 

75 

-120 

126 

25 

-167 

19 

1800 

83 

-100 

134 

43 

-171 

25  AUG 

0700 

67 

-123 

135 

61 

-156 

21 

1800 

SO 

-85 

139 

93 

-142 

24 

26  AUG 

0800 

63 

-109 

140 

100 

-131 

20 

1800 

75 

-81 

144 

114 

-129 

23 

27  AUG 

0700 

65 

-92 

150 

120 

-123 

22 

1900 

66 

-65 

149 

127 

-109 

24 

28  AUG 

0800 

60 

-57 

150 

131 

-98 

20 

. 1 "  e0600 

1800 

71 

-31 

154 

141 

-103 

19 

29  AUG 

0800 

59 

-16 

157 

143 

-97 

13 

1 800 

79 

18 

160 

153 

-78 

17 

30  AUG 

0800 

64 

22 

159 

150 

-87 

13 

1800 

79 

39 

160 

157 

-73 

17 

31  AUG 

0800 

67 

36 

159 

153 

-76 

15 

2000 

75 

45 

159 

157 

-66 

18 

1  SEP 

0800 

68 

44 

159 

155 

-77 

17 

.  2  " 

1200 

68 

43 

157 

154 

-74 

18 

1 600 

69 

44 

156 

153 

-60 

18 

1800 

70 

44 

155 

153 

-67 

17 

.  1  '■ 

2  SEP 

0800 

61 

38 

144 

147 

-167 

17 

.  4" 

1 800 

78 

1 56 

155 

-202 

18. 

Table  12 

SP  Values  in  Millivolts  for  Long-Term  Monltorln] 
Array  of  Electrodes 


ELECTRODES  SOIL  TEMP 


DATE 

TIME 

#1 

#2 

#3 

#4 

#5 

6" 

3  SEP 

0800 

108 

34 

180 

176 

-205 

18 

1800 

130 

26 

197 

191 

-204 

18 

4  SEP 

0700 

1 33 

11 

202 

197 

-208 

18 

1800 

146 

0 

213 

206 

-208 

19 

5  SEP 

0800 

149 

-49 

217 

210 

—212 

18 

2100 

164 

-78 

231 

219 

-178 

19 

6  SEP 

0900 

167 

-87 

238 

228 

-198 

18 

2000 

172 

-77 

238 

229 

-182 

19 

7  SEP 

0800 

158 

-85 

230 

228 

-179 

17 

1800 

161 

-65 

228 

227 

-150 

19 

8  SEP 

0800 

157 

-30 

230 

232 

-133 

18 

1800 

164 

—3 

230 

235 

-114 

18 

9  SEP 

0800 

146 

-26 

223 

229 

-105 

16 

1900 

147 

1 

217 

225 

-101 

17 

10  SEP 

0700 

143 

1 

'70'? 

229 

-92 

16 

1800 

142 

27 

220 

229 

-76 

18 

11  SEP 

0800 

129 

21 

“7 '7 '7 

230 

-76 

18 

1800 

134 

39 

221 

229 

-71 

19 

12  SEP 

0800 

126 

44 

222 

230 

-75 

19 

1800 

131 

57 

223 

231 

-76 

19 

13  SEP 

0900 

119 

63 

222 

230 

-64 

18 

1800 

127 

76 

220 

228 

-51 

19 

14  SEP 

0800 

115 

78 

221 

228 

-61 

16 

1700 

127 

88 

221 

230 

-45 

18 

15  SEP 

0700 

113 

82 

216 

224 

-61 

17 

1800 

119 

89 

217 

226 

-52 

19 

16  SEP 

0800 

114 

86 

217 

225 

-59 

18 

1800 

119 

91 

216 

225 

-41 

19 

17  SEP 

0700 

110 

87 

217 

225 

-60 

15 

18  SEP 

0800 

118 

91 

221 

231 

-51 

15 

1800 

126 

97 

222 

232 

-34 

17 

19  SEP 

0800 

120 

92 

221 

229 

-50 

17 

1900 

123 

92 

218 

227 

-44 

18 

20  SEP 

0900 

123 

93 

222 

230 

-52 

17 

21  SEP 

0800 

121 

92 

221 

230 

-61 

17 

22  SEP 

0700 

127 

98 

225 

234 

-59 

18 

23  SEP 

0800 

131 

100 

228 

239 

-62 

18 

2000 

139 

106 

234 

234 

-54 

20 

24  SEP 

0700 

139 

108 

239 

238 

-58 

19 

1900 

121 

87 

217 

215 

-76 

20 

25  SEP 

1900 

126 

89 

218 

216 

-76 

21 

26  SEP 

2000 

127 

88 

218 

217 

-78 

21 

27  SEP 

1400 

129 

90 

221 

222 

-80 

21 

28  SEP 

1  400 

134 

90 

225 

225 

-80 

21 

30  SEP 

0800 

126 

88 

226 

231 

-114 

19 

1  OCT 

1100 

141 

97 

236 

242 

-105 

21 

2  OCT 

1700 

145 

94 

233 

239 

-109 

21 

3  OCT 

1200 

144 

94 

236 

243 

-112 

19 

6  OCT 

0800 

139 

87 

235 

242 

-128 

18 

7  OCT 

0800 

134 

85 

233 

243 

-129 

16 

a  OCT 

0900 

132 

88 

234 

247 

-129 

IS 

1 800 

147 

95 

241 

253 

-110 

18 

RAINFALL 


tl 


II 


Table  13 

Record  of  Precipitation  from  Rain  Gage  at  Long-Term 
SP  Electrode  Monitoring  Station 


DATE 

INCHES 

7-21-86 

.  5 

7-22-86 

.  1 

7-23-86 

.  4 

7-26-86 

.25 

8-06-86 

.8 

8-17-86 

.  5 

8-19-86 

.6 

8-20-86 

1.  1 

8-21-86 

.2 

8-28-86 

.  1 

9-01-86 

9-02-86 

•  5 

9-03-86 

9-04-86 

.  1 

9-0S-86 

10-13-86 

.25 

10-14-86 

.25 

10-25-86 

1 . 35 

11-01-86 

.5 

11-05-86 

1.2 

11-07-86 

.4 

11-11-86 

.  1 

11-18-86 

.25 

11-20-86 

.32 

11-26-86 

.8 

12-02-86 

1.5 

12-09-86 

.  3 

12-1 1-86 

.9 

12-18-86 

.4 

12-24-86 

2.4 

01-01-87 

.4 

01-15-87 

.  15 

01-19-87 

1.5 

01-30-87 

r 

01-31-87 

1 

02-12-87 

.  1 

02-28-87 

n 

03-08-87 

.45 

03-16-87 

.5 

03-19-87 

.4 

03-27-87 

.  3 

03-30-87 

1.2 

03-31-87 

.6 

04-03-87 

.  5 

04-15-87 

2.5 

04-16-87 

2.9 

04-17-87 

1.4 

04-24-87 

3.  6 

05-04-07 

1.  1 

NOTES 


as  snow 


snow  melt 
snow  melt 


as  snow 


as  snow 


Table  14 

Record  of  SP,  Soil  Temperature,  Soil  Moisture,  and  Standard  Deviation  at  Long-Term 
Monitoring  Site  for  October  to  November  1986  (Data, reduced  to  24-hr 

average  from  ADC  record) 
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Table  15 

Record  of  SP,  Soil  Temperature,  Soli  Moisture,  and  Standard  Deviation  at  Long-Term 
Monitoring  Site  for  December  1986  (Data  reduced  to  24-hr  average  from  ADC  record) 
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Table  16 

Record  of  SP,  Soil  Temperature.  Soil  Moisture,  and  Standard  Deviation  at  Long-Term 
Monitoring  Site  for  January  1987  (Data  reduced  to  24-hr  average  from  ADC  record) 
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Record  of  SP,  Soil  Temperature,  Soil  Moisture,  and  Standard  Deviation 
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Record  of  SP,  Soil  Temperature,  Soil  Moisture,  and  Standard  Deviation  at  Long-Term 
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Table  20 

SP  Values  (millivolts)  at  the  Moore  Sinkhole  A,  July  to  October  1986 
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SP  Values  (millivolts)  at  the  Three  Moore  Sinkholes  from  January  to  April  1987 
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Figure  1.  Schematic  of  flow- induced  negative  e-treaming  potentials  (from 

Cooper  1983) 
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Figure  2.  Topographic  map  of  study  area  (contours  in  metres). 
A  =  McHenry  sinkhole  A;  B  =  McHenry  sinkhole  B;  C  =  Whose  Hole 
Cave;  D  *  Electrode  monitoring  site;  E  =  Spring  Branch 
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Figure  3.  Polarization  equilibrium  of  electrodes  following 
initial  placement  at  the  electrode  monitoring  test  site 
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Figure  6.  SP  for  electrode  array  for  two  different  reference 
electrodes,  3  August  1986,  Spring  Branch  site 
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Figure  8.  SP,  streamflow,  and  rainfall.  Spring  Branch  site.  SP  data 
show  different  magnitudes  for  the  two  halves  of  the  array,  relative 
to  their  position  in  relation  to  the  stream  sink  (swallet) 
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Figure  10.  Map  with  SP  contours  (In  millivolts)  and  electrode 

array.  Ford  site 
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SP  for  an  anomalous  electrode  (No.  9),  streamflow,  and  temperature.  Ford  site 


Figure  12.  Topographic  map  of  study  area  (contours  In  metres)  show¬ 
ing  relationship  between:  A  (McHenry  sinkhole  A),  B  (McHenry  sink¬ 
hole  B) ,  C  (Whose  Hole  Cave),  and  D  (electrode  monitoring  site) 
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Figure  20.  Array  of  12  electrodes  at  entrance  sinkhole  over  Whose 
Hole  Cave  (Arrows  indicate  direction  of  intermittent  seepage  flow 

into  and  through  the  cave) 
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Figure  29.  ACD  data  from  electrode  monitoring  site  showing  soil 
temperature,  soil  moisture,  and  SP  response  of  electrode  No.  3 
(in  streambed)  to  rainfall  of  2.5  in.  in  24  hr 
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Figure  31.  Change  In  SP  Induced  by  water  discharge  tests  at 

Moore  sinkhole 


Figure  34.  SP  at  Moore  sinkhole  showing  change  in  SP  induced  by  water 
discharge  test  and  the  influence  of  temperature,  18  and  19  October  1986 
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TO-metre  survey 

Figure  36.  EM  conductivity  surveys  at  Moore  sinkhole  A,  December  1986, 
contours  in  nunhos/m  (Circles  are  SP  electrode  arrays) 


Figure  37.  EM  conductivity  surveys  at  Moore  sinkhole  B,  December  1986, 
contours  in  nanhos/m  (Circles  represent  SP  electrode  arrays) 


Figure  38.  EM  conductlvl 
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Figure  39.  Moore  site  sinkholes  A,  B,  and  C.  Dashed  lines  are 
topographic  contours;  solid  lines  are  conductivity  contours  In 
mmhos/in  from  10  m  survey  in  December  1986  (Circles  represent  the 
diameter  of  outer  ring  of  SP  electrodes  Installed  In  each 

sinkhole) 


Figure  42.  SP  contours  In  millivolts  at  Moore  sinkhole  B, 
January  1987  (Circles  represent  electrode  arrays) 
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Figure  44.  SP  contours  in  millivolts  at  Moore  sinkhole  A,  SP  values 
averaged  for  period  July  -  October  1986  (Circles  represent 

electrode  arrays) 
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Figure  45.  Dally  average  SP  for  electrode  No.  1  at  long-term  monitoring  site  over  7-month  period 
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Figure  48.  Daily  average  SP  for  electrode  No.  4  at  long-term  monitoring  site  over  7-month  period 


Figure  49.  Daily  average  SP  f< 


Figure  50.  Soil  temperature  (degrees  C)  at  depth  of  20-cm  at  the  long-term  monitoring  site 
for  a  period  of  7  months.  Data  are  daily  averages  from  ADC  record 
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Figure  51.  Dally  average  soli  moisture  at  20-cm  depth  at  long-term  monitoring  site  over  7-month  period 
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Figure  52.  Relationship  of  rainfall  to  SP  at  electrode  No.  3  (in  bed  of  the  intermittent  stream) 
at  long-term  monitoring  site.  Sharp  drops  in  SP  occur  when  the  stream  was  flowing 


53.  Correspondence  of  soil  temperature  and  inverted  SP  record  of  electrode  No. 

at  the  long-term  monitoring  site 
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